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Abstract. We address the degree-diameter problem for Cayley graphs of Abehan groups 
(Abehan graphs), both directed and undirected. The problem turns out to be closely related 
to the problem of finding efficient lattice coverings of Euclidean space by shapes such as 
octahedra and tetrahedra; we exploit this relationship in both directions. For 2 generators 
(dimensions) these methods yield optimal Abelian graphs with a given diameter k. (The 
results in two dimensions are not new; they are given in the literature of distributed loop 
networks.) We find an undirected Abelian graph with 3 generators and a given diameter k 
which we conjecture to be as large as possible; for the directed case, we obtain partial results. 
These results are connected to efficient lattice coverings of by octahedra or by tetrahedra; 
computations on Cayley graphs lead us to such lattice coverings which we conjecture to be 
optimal. (The problem of finding such optimal coverings can be reduced to a finite number of 
nonlinear optimization problems.) We discuss the asymptotic behavior of the Abelian degree- 
diameter problem for large numbers of generators. The graphs obtained here are substantially 
better than traditional toroidal meshes, but, in the simpler undirected cases, retain certain 
desirable features such as good routing algorithms, easy constructibility, and the ability to 
host mesh-connected numerical algorithms without any increase in communication times. 



Introduction. The degree-diameter problem for graphs is the following question: What 
is the largest number of vertices a graph (undirected or directed) can have if one is given 
upper bounds on the degree of each vertex and on the diameter of the graph (the maximum 
path-distance from any vertex to any other)? One application for such graphs is in the 
design of interconnection networks for parallel processors, where one wants to have a large 
number of processors without requiring a large number of wires at a single processor or 
incurring long delays in messages from one processor to another. For more information on 
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the (undirected) degree-diameter problem, see Dinneen and Hafner [10]; up-to-date results 
can be found at the following Web site: 

http : //www-mat .upc . es/grup_de_graf s/graf s/taula_delta_d.htinl 

A desirable extra property of such networks is that they look identical from any pro- 
cessor; this means that the graphs one uses should be vertex-transitive, i.e., for any two 
vertices x and y there is an automorphism of the graph which maps x to y. Here we will 
restrict ourselves to a special class of vertex-transitive graphs, called Cayley graphs. A 
Cayley graph is specified by a group and a set of generators for this group; the vertices 
of the graph are the elements of the group, and there is an edge from a: to y if and only 
if there is a generator g such that y = xg. (It can be shown that every vertex-transitive 
graph is isomorphic to a generalized form of Cayley graph called a Cayley coset graph [19]. 
In this paper, though, we will look only at Cayley graphs. For Abelian groups, this is 
no loss of generality since every Cayley coset graph of an Abelian group is isomorphic 
to a Cayley graph of an Abelian group.) In a directed Cayley graph from a group on d 
generators, every vertex has in-degree and out-degree d; if d generators are used to form 
an undirected Cayley graph, then the degree of each vertex is the number of generators 
of order 2 plus twice the number of generators of order greater than 2 (unless there are 
redundant generators). So we will usually talk about Cayley graphs on a given number 
of generators rather than of a given degree; the cases where some generators have order 2 
and hence only contribute 1 rather than 2 to the degree of an undirected Cayley graph 
will be handled separately. 

The most straightforward approach for trying to find large Cayley graphs of small 
diameter on a given number d of generators is to examine various groups, look at some or 
all possible sets of d generators for such a group, and check whether each such set in fact 
generates the group and, if so, what the diameter of the graph is. But this can be a very 
large task even for relatively small groups and generating sets. In this paper, we will use 
a different approach which facilitates studying many groups and generating sets at once; 
it yields provably optimal results for some families of groups, and good lower and upper 
bounds for others. 

In its most general form, the idea is as follows. Let F be the free (universal) group 
on d generators. Then, for any group G and any set of d generators for G, there is 
a homomorphism -k: F G which maps the canonical generators for F to the given 
generators for G; clearly tt is surjective. Let be the kernel of tt. Then is a normal 
subgroup of F, and \F : N\ = \G\] in fact, G is isomorphic to F/N , and the Cayley graph 
of G with the given generators is isomorphic to the Cayley graph of F/N with the canonical 
generators for F. Let S be the set of elements of F which can be expressed as a word of 
length at most k in the generators. (If one is interested in undirected Cayley graphs, then 
such words may use inverse generators g~^ as well as generators; for the directed case, 
only words using generators and not inverse generators are allowed.) 

Proposition 1. The Cayley graph for G on the given generators has diameter at most k 
if and only if SN = F . 
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Proof. First, suppose SN = F. Let a be an arbitrary element of G; then a = n{x) 
for some x, and x can be written in the form wy with w & S and y & N. Hence, 
a = iriwy) = 'n{w)'K{y) = 7r{w). Now w can be written as a word of length at most k in 
the generators of F, so a = 7t{w) can be written as the same word in the corresponding 
generators of G. Since a was arbitrary, the Cayley graph has diameter at most k. 

Now suppose the Cayley graph has diameter at most k. Let x be any element of F; 
then 7r{x) can be written as a word w' of length at most k in the generators of G. Let w 
be the corresponding word in the generators of F; then 7t{w) = w' = 7r{x), so 7r{w~^x) is 
the identity of G, so w~^x e N. Hence, x = w{w~^x) e SN, as desired. ■ 

So finding a Cayley graph on d generators with diameter k whose size is as large as 
possible is equivalent to finding a normal subgroup N of F such that SN = F and | F : | 
is as large as possible. Of course, we immediately get the upper bound |F : < l^l, but 
this is probably not attainable. 

Unfortunately, the collection of normal subgroups of F is so large and varied as to be 
unmanageable. So what we will do instead is restrict ourselves to certain families (usually 
varieties) of groups; this allows us to replace F with a free group for the family in question, 
which may be much easier to work with. For instance, if we only consider the Cayley graphs 
of Abelian groups, then we can replace F with the free Abelian group on d generators, 
which is simply Z^; the normal subgroups of are well understood and relatively easy to 
work with. It turns out that this reduces the degree-diamctcr problem for Abelian Cayley 
graphs to interesting problems about lattice coverings of Euclidean space by various shapes. 
Some of these problems can be solved completely, giving optimal Abelian Cayley graphs; 
others are still open. 

A simple path-counting argument gives upper bounds for the size n of a Cayley graph 
with d generators and diameter limit k: in the directed case, 

n<l + d + d'^-\ l-d'' = — 



and in the undirected case, 

n<l + 2d + 2d{2d -!) + ••• + 2d{2d - 1)'=-^ = "^^^"^ ~ ^ . 

Ct J- 



(The formulas for d = 1 are A; + 1 and 2k + 1.) These limits arc well-known and actually 
apply to the degree- diameter problem for arbitrary graphs; they are known as the Moore 
bounds. For d = 1, these limits are attained by simple cycle graphs, which are Cayley 
graphs of cyclic groups. 

In most cases, we will find that the attainable values for n using Cayley graphs in the 
families we consider do not approach these upper bounds; the equations defining the fami- 
lies force many paths to be redundant. But the extra structure provided by the groups may 
provide compensating advantages in parallel computers, such as good routing algorithms, 
easy constructibility, and the ability to map common problems onto the architecture. In 
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particular, many of the Cayley graphs of Abehan groups that we discuss in this paper are 
multi-dimensional rectangular meshes with additional connections at the boundary. Thus, 
mesh calculations with natural boundary conditions are trivially mapped into these graphs, 
while the extra connections are utilized only when global communications are being carried 
out. In addition, the mesh nature of these graphs allows the physical construction of the 
network to be carried out with relatively short wires. This will be discussed further below. 

In a separate paper, we will examine some other varieties of groups for which similar 
analyses of Cayley graphs are feasible. 

We would like to (and hereby do) thank Michael Dinneen, whose computer work deter- 
mining the diameters of many Cayley graphs yielded some of the specific groups listed here 
and suggested promising families of groups to examine. Thanks also to Francesc Comellas, 
for pointing us to the existing literature on directed loop networks (in particular, the very 
useful survey paper of Bermond, Comellas, and Hsu [3]), and to Alexander Hulpke for 
discussions on automorphisms of Abelian groups. 

Abelian groups. In the rest of this paper, we will examine the Cayley graphs arising 
from Abelian groups. Toroidal meshes and hypercubes are examples of such graphs. 

The degree- diameter problem for Abelian Cayley graphs has been considered by others. 
In particular, Annexstein and Baumslag [2] show that the number of generators d, diameter 
k, and size n of a directed Abelian Cayley graph satisfy 

k > n{n^/'^); 

in fact, a d < n^/*^, then 

k > n{dn^/'^). 

They also discuss similar results for Cayley graphs of nilpotent groups. 

In addition, Chung [6] has constructed directed Abelian Cayley graphs G with n = p* — 1, 
d = p, and 

r mogt 

k< 2t+- — ^- - 

logp-21og(t-l) 

for any positive integer t < ^/p + 1, where p is a prime. (Note that Chung's examples 
have diameters which are small compared to the number of generators; in contrast, we will 
concentrate here on graphs with a small fixed number of generators and relatively large 
diameters.) Her methods involve estimates of the second eigenvalue of the Laplacian of the 
adjacency matrix of G. For more about estimating the diameter of general graphs from 
knowledge of this eigenvalue, see Chung, Faber, and Manteuffel [7]. This eigenvalue is also 
connected to the sphere packing problem for real lattices; see Urakawa [21]. 

The more specific case of Cayley graphs of cyclic groups has been studied more exten- 
sively; such graphs are usually referred to by some variant of the phrase 'loop networks.' 
The survey paper of Bernard, Comellas, and Hsu [3] is an excellent guide to the literature 
in this area. 

We start by taking care of some generalities and notational matters. In this paper we 
will use the symbol -|- for the group operation(s). Let be the cyclic group of order m 
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(for definiteness, the set {0, 1, . . . , m — 1} with the operatfon of additfon modufo m). The 
groups and the infinite cychc group Z each have a canonical generator, 1; but they 
also have other sets of generators, and some of these will be important later. 

When one has groups Gi, G2, ■ • • , and a set of generators for each, then one can get a 
set of generators for the product Gi x G2 x ■ ■ ■ x by putting together the given generator 
sets. More precisely, for each i < I and each generator g of Gi, let ei{g) be the element of 
the product group which has the identity element of Gj as its j'th coordinate for all j < I 
except i; the i'th coordinate is g. Then the set of all elements ei{g) is a natural generating 
set for the product group. (The resulting diameter for the product group is the sum of the 
diameters of the groups Gj.) In the case where the groups G, are cyclic groups with the 
canonical single generators, we write simply for the Z-tuple with 1 at the i'th coordinate 
and elsewhere. 

A two-dimensional toroidal mesh is simply the Cayley graph of the group x Z„ 
with the canonical generators ei = (1,0) and 62 = (0, 1); higher-dimensional meshes are 
obtained from longer products. For the two-dimensional case, the number of vertices is mn, 
the degree is 2 in the directed case and 4 in the undirected case (assuming m, n > 3), and 
the diameter is m + n — 2 in the directed case, [m/2j -|- [n/2j in the undirected case. The 
calculations in three or more dimensions arc analogous. 

The (i-dimensional hypercube is the Cayley graph of the group Z2 with the canonical 
generators; since one gets bidirectional links in any case, we may as well just talk about 
the undirected version. In this case the size of the graph is exponential in the diameter, 
but only because the degree also grows with d: the size is 2^^, the degree is d (not 2d, 
because the generators have order 2), and the diameter is also d. 

We will see that with a fixed small number d of generators, one can obtain nearly optimal 
results for undirected Abelian Cayley graphs by using a twisted toroidal mesh; the twist 
allows one to multiply the number of nodes in the ordinary toroidal mesh by 2^~^ without 
increasing the diameter. For two dimensions we can get exactly optimal results by slightly 
adjusting this graph; for higher dimensions, the optimal size is still open. For the directed 
case, we can again get optimal results in two dimensions, but the higher-dimensional case 
is again unsolved (and rather strange even in three dimensions). 

To get these results, we argue as in Proposition 1, but specialize to the case of Abelian 
groups. Hence, instead of the free group on d generators, we can use the free Abelian 
group on d generators, which is simply 7j!^ with the canonical generators e^, 1 < i < For 
any Abelian group G generated by gi, . . . there is a unique homomorphism from Z^ 
onto G which sends to gi for all i. Let be the kernel of this homomorphism; then G 
is isomorphic to Zi^/N, and the Cayley graph of G with the given generators is isomorphic 
to the Cayley graph of Z^/N with the canonical generators for Z^. 

Given a diameter limit /c, let Sk be the set of elements of Z^ which can be expressed as 
a word of length at most k in the generators e^, which are allowed to occur positively or 
negatively. (The dimension d will be clear from context.) Then can also be described 
as the set of points in Z*^ at distance at most k from the origin under the (Manhattan) 
metric: 

Sk = {{xi,.. .,Xd) e Z^: |a;i| + \- \xd\ < k}. 
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Let S'^. be the subset of Sk consisting of those elements whose coordinates are all nonneg- 
ative; these are the elements which can be expressed as words of length at most k in the 
generators where only positive occurrences of the generators are allowed. Then Sk looks 
like a regular dual d-cube (a square for d = 2, an octahedron for d = 3), while S'^ looks 
like a right d-simplex (a triangle for (i = 2, a tetrahedron for d — 2>). 
Now, by the same proof as for Proposition 1, we get: 

Proposition 2. Let G, N , and gi, . . . , be as above. Then the undirected Cayley graph 
for G and gi, . . . ,g(i has diameter at most k if and only if Sk + N = Zi^, and the directed 
Cayley graph for G and gi, . . . ,g(i has diameter at most k if and only if S'j. -\- N = 7/^ . ■ 

So \Sk\ and jS'^j give upper bounds for the undirected and directed versions of this case 
of the degree-diameter problem. It is not hard to show that 

,.;,=er), 

so \S'f^\ = k'^/dl + 0{k'^~^) for fixed d. For \Sk\, we easily get the asymptotic form I^Sfcl = 
/d\ -\-0{k'^~^) for fixed d, but exact formulas are harder; Stanton and Cowan [20] give 
several, such as 

In particular, when d is 1, 2, or 3, the formula for \Sk\ is 2k + 1, 2k^ + 2A; + 1, or (4A;^ + 
+ 8/c + 3)/3, respectively. 

Lattice coverings and tilings. Proposition 2 tells us that, to find an optimal undirected 
(directed) Cayley graph of diameter k on d generators, we should look for a subgroup N 
of such that Sk + N (S*^ + N) is all of and the index [Z*^ : A?"! is as large as possible; 
the largest index we can hope for is |5'fc| (|5'^|). But the structure of subgroups of Z'^ 
of finite index (which are all normal, of course) is well known; they are precisely the d- 
dimensional lattices in Z*^. Because of this, we will use the letter L instead of N for such 
subgroups of Z'^ for the rest of this paper. 

A d-dimensional lattice L in Z^ is specified by d linearly independent vectors vi , . . . , 
in Z'^; L is the set of all integral linear combinations of these vectors. We have \Z : L\ = 
I det M|, where M is the d x d matrix whose z'th row is v^, for z = 1, . . . , d. 

Note that any bounded set contains only finitely many members of L. It follows that, 
if 5" is a bounded subset of Z^ and x is a point in Z^, then there are only finitely many 
V e L such that x e S + v. 

Also note that L, or indeed the entire group Z^, has a linear ordering -< which is 
compatible with addition: x -< y implies x + z -< y + z. To define first choose a 
direction (a nonzero vector v in R'^), and put x -< y if y is farther in this direction 
than x is (x ■ V < y ■ v). If two vectors are at the same distance in this direction, then 
compare them in a second direction; repeat until all ties are broken. One example of this 
is lexicographic order: compare according to the first coordinate, then according to the 



THE DEGREE-DIAMETER PROBLEM FOR ABELIAN CAYLEY GRAPHS 



7 



second coordinate if the first coordinates are equal, and so on. Or, in this discrete case, 
one can choose the initial direction so that it distinguishes all points and no tie-breaking 
is necessary; for instance, if v = (1, tt, tt^, . . . , tt'^"^), then we never have x • v = y • v for 
distinct x, y in Z'^, so we can just define x ^ y to mean x ■ v < y ■ v. 

A lattice covering of Z*^ by a set 5" C Z*^ is a collection of translates of S by members 
of a lattice L (i.e., {5 + v: v e L}) which covers Z^. If the translates are disjoint, so that 
each point of Z*^ is covered exactly once, then we have a lattice tiling of Z*^ by S. 

If we have a lattice covering as above, then j^l > \Z^ : L|; if it is a tiling, then l^l = 
IZ'^ : L|. So we can measure the extent to which a covering is 'almost' a tiling by one of 
two numbers: the density of the covering, which is jZ'^ : L| /l^"! > 1 (this is the average 
number of sets in the covering to which a random point of Z^ belongs), or the efficiency 
of the covering, which is l^l/ jZ'' : L| < 1. So Proposition 2 tells us that, in order to get 
the best possible Abelian Cayley graph on d generators with diameter k, we must find a 
lattice covering of Z'^ by Sk or S'^. whose density (efficiency) is as small (large) as possible. 
We now give one more reformulation of the question. 

Lemma 3. Suppose we have a lattice covering of Z*^ using a bounded set S and the 
lattice L. Then there is a set T C S such that the translates of T by L form a lattice tiling 
of Zd. 

Proof. Let -< be a linear order of L compatible with addition. Now let T be the set of all 
points in S which are not in any of the sets S' + vforveL,v:^0. We will see that every 
point X is in exactly one of the sets T + v for v e L. 

Fix X. As noted before, x is in only finitely many of the translates S' + v, so let w be 
the -<-greatest member of L such that x G 5" + w. Let y = x — w G 5". Then y cannot be 
in + V for V :^ in L, because, if it were, we would have x^y + wGS'-l-v-l-w and 
V + w :^ w, contradicting the maximality of w. So y G T and x G T + w. 

Now suppose x = y-|-w = y' + w' where w and w' are distinct members of L and 
y and y' are in T (and hence in S). Then v = y — y' = w' — w is a nonzero member of L, 
so either v or v -< 0. In the former case, y = y' + vG5'-|-v contradicts y G T; in the 
latter case, y'=y — vg5' + (— v) contradicts y' G T. ■ 

Note that such a set T must be of cardinality \Z'^ : L\, which is the size of the Cayley 
graph of Z'^/L. So the size of the largest undirected (directed) Abelian Cayley graph on 
d generators with diameter k is equal to the size of the largest subset T of {S^,) such 
that there is a lattice tiling of Z^ using T. 

Approximation by lattice coverings of real space. The study of lattice coverings 
and lattice tilings is more familiar for R'' than for Z^. We will show that real coverings 
can be approximated to some extent by integer coverings, and vice versa, so that known 
results from the real context can be transferred to the integer lattices we are interested in. 

The definitions of lattice, lattice covering, and lattice tiling are the same in R'^ as in Z^, 
except that we allow boundaries to be shared in the definition of a tiling. This lets us work 
throughout with closed sets (usually polyhedra with their interiors) instead of having to 
keep some of the boundary points and discard others. Most of the results above for integer 
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lattices go through verbatim for real lattices, including Lemma 3 and its proof (although 
we will probably use the closure of T rather than T itself in practice). The main difference 
is that the absolute determinant | det M \ of the matrix formed from the generators of a 
lattice L is not IR'^ : L\ (which is infinite). Instead, this determinant is the measure of 
the parallelepiped determined by the generating vectors; this parallelepiped gives a lattice 
tiling of using the lattice L. It follows easily that any other set S which gives a lattice 
tiling using the lattice L must have measure |detM| (barring pathological cases of non- 
measurable sets or positive- measure boundaries). Such a set S is called a fundamental 
region for the lattice L. One can now define the density or efficiency of a covering by 
dividing the measure of the covering set by this determinant or vice versa. 

One can transform a lattice covering of 7/^ using S into a lattice covering of R*^ by 
replacing each point of S with a unit cube (i.e., replace S with S + U where U is a fixed 
unit d-cnhe with edges parallel to the coordinate axes); the two coverings will have the 
same density. However, transforming results in the other direction is harder, because 
the real results usually involve actual triangles, octahedra, etc. rather than polycube 
approximations, and the lattices used will often not be integer lattices. We will now 
present results that allow us to get around these difficulties. 

In R'', let Sk be the closed Z^-ball of radius k at the origin: 

Sk = {{xi, Xd): \xi\-\ 1- \xd\ < k]. 

Let S'^ be the set of nonnegative points in Sk'- 

S'k = {{xi, ...,Xd):xi,...,Xd>0,xi-\ \-Xd<k}. 

Let L be any lattice in R'^. 

Proposition 4. 

(a) If Sk + L covers TJ^ , then S^j^dii + covers R'^. 

(b) If S'^+L covers TJ^ , then S'-^^^ + L covers R'^. 

Proof, (a) By the triangle inequality for distance, we have Sk + Sd/2 ^ Sk+d/2: so 
Sk + L + Sd/2 Q Sk+d/2 + L; therefore, it suffices to show that 7/^ + Sd/2 = For any 
X e R*^, let y be the element of Z*^ nearest to x (i.e., round each coordinate of x to the 
nearest integer); then x — y is in S'd/2, so x is in Z*^ + Sd/2- 

(b) Similarly, this follows from the fact that 7/^ + 3'^, = R*^, which is proved in the same 
way as above (round each coordinate of x downward instead of to the nearest integer) . ■ 

One can argue in the same way within 7/^ to get: 

Proposition 5. If L is a lattice in 7'^ and m is a positive integer, then: 

(a) If Sk + L covers 7^, then Smk+lm/2id + itT'L covers 7^. 

(b) If S'f. + L covers 7^, then S'^kMm-\\d + "^-^ covers 7'^. 
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Proof. For (a), clearly Smk+nT'L covers mZ'^, so, as in the preceding proposition, it suffices 
to note that mZ^ + Sim/2]d covers Z*^, because we can just round any member of Z^ to 
the nearest member of mZ'^. Similarly, (b) holds because mZ^ + S'^^_^s^^ covers Z*^. ■ 

Of course, one can get similar results for sets other than Sk and S'f^. 

Using Proposition 4, it is easy to get from a covering of Z*^ using an integer lattice to a 
covering of R'' using a real lattice (and using the real shape 5"^ or 5"^); if k is large relative 
to d, then the two coverings have about the same efficiency. We will now show that one 
can move in the other direction as well. 

First, we give a useful criterion for deciding whether one has a lattice covering of R'^. 

Proposition 6. Suppose S is a nonempty subset of R^ and L is a lattice in R'^. If there 
is an e > such that S + L covers all points within distance s of S, then S + L covers R''. 

Proof. Clearly there is some point xq in S' + L. We will show that, if x e S+L and 
the distance 5(x, y) is less than e, then y ^ S + L. Applying this once shows that all 
points within distance e of xq are in 5" + L; applying it again shows that all points within 
distance 2e of xq are in S + L; since this can be repeated forever, we find that all points 
of R'^ are in S + L. 

Let X and y be as above. Find v e L such that x e S + v. Then y — v is within 
distance £ of x — v G S', so there is v' e L such that y — v e S + v'. Hence, y e 5' + v' + v, 
so y e 5' + L, as desired. ■ 

Using this, we can now show that, if one has a lattice covering using a bounded subset 
of R^, then one can perturb the lattice slightly and still get a lattice covering using a 
slightly larger subset of R''. 

Proposition 7. Let S be a bounded subset of H'^, and let L be a lattice in H'^ such that 
S + L = H'^; let Vi, . . . , Vrf be a list of generators for L. Then there are positive numbers 

1] and p such that, for all r G (0, 1), if the distance (5(vi, vQ is less than rrj for all i < d, 
then S~^ + L' = H'^ , where L' is the lattice generated by v^, . . . , and S^ is the set of 
points within distance rp of S. 

Proof. The number of members of L within any bounded part of R^ is finite, so it only 
takes finitely many translates of S by members of L to cover any bounded part of R'^. 
In particular, there is a number M > such that the sets S + aiVi + • • • + a^Vd for 
(ai, . . . , ad) G Z^ with |ai| + • • • + \ad\ < M cover all points within distance p of S, for 
some p > 0. Let rj = p/M. 

Let r be any positive number less than 1; we must see that, if S~^ and L' are defined 
as above, then 5+ + L' — R^. By the preceding proposition, it will suffice to show that 
S'^ + L' covers all points within distance (1 — r)p of 5'+. Suppose y is within distance 
(1 — r)p of S~^] then y is within distance (1 — r)p + rp — p of S, so there exist integers 

oi, . . . , Od with |ai| H h \ad\ < M and a point x G 5 such that y = x + aiVi H l-ad'^d- 

Let x' = X + Yli=i - vO; then 

d d 

5(x,x ) < ^ \ai\6{vi,v[) < ^ \ai\rr] < Mrr] = rp. 
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SO x' e S'^. Since y = x' + aiv'j^ + h a^v^, we have y e 5'+ + L', as desired. ■ 

Note that, in the above propositions, 'distance' need not be Euclidean distance; it can 
be any metric arising from a norm on R''. For our present purposes, it will be most 
convenient (but not essential) to use Z°°-distance: 5{x,y) = maxj \xi — yi\. 

Theorem 8. Suppose one has a lattice L in R*^ such that Sk + L covers R'^; let Vi, . . . ,Vd 

be generators for L. Then there is a constant c such that, for all sufficiently large real 
numbers t, if is obtained from tv^ by rounding all coordinates to the nearest integer, 
and L is the lattice generated 6y wi, . . . , w^, then Stk+c+L covers IV^. The same statement 
holds for S[. and S[f^j^^ instead of Sk and Stk+c- 

Proof. For the Sk case, let S = Sk and find i] and p as in the preceding proposition, letting 
the distance 6 be the l°° metric. Let c be any fixed number greater than dp/2r]. Then, 
for any t > c/dp, if we let r = c/dpt, then r < 1 and l/2t < rij. If we define Wj and L as 
above, and let v'- = Wi/t and L' = L/t, then (5(vi, vQ < l/2t for each i, so we can conclude 
that S~^ + V covers R'^, where S^ is the set of points within distance rp of Sk- It is easy 
to see that S^ C Sk+drp, so Sk+drp + L' covers R^. Hence, Stk+tdrp + tL' covers R'^; but 
tdrp — c and tL' — L, so we are done. 

The proof for S'^. is almost the same. Let S = S'^ and apply the preceding proposition 
(using the l°° metric again) to get rj and p. Fix c > dp/rj. For any t > c/2dp, if we let 
r = c/2dpt, then r < 1 and l/2t < rrj. Now define w^, L, v^, L', and 5'"'" as above, and 
conclude again that 5"+ + L' covers R^. One can easily check that 5"+ C S'j^j^2drp ~ ^P^' 
where u = (1,...,1). Hence, S'^j^2drp " + ^' covers R"^, so 5'^+2drp + ^' covers 
pj^d _|_ _ j^d_ ]\Jq^ multiply by t to see that S[j^j^^ + L covers R''. ■ 

Again, it is easy to modify this proof to work for other sets in place of Sk or S'j^. Also, 
the proof is quite effective, allowing one to compute specific values of c and t which work 
for a given lattice L (assuming it is feasible to compute M and p). 

The covering Stk+tL has the same efficiency as the covering Sk + L; since Stk+c + L is a 
relatively slight perturbation of Stk + tL when t is large, it has almost the same efficiency. 
Since L is an integer lattice, the fact that Stk+c + L covers R*^ implies that S^tk+c} + L 
covers Z^; again the efficiency is almost the same if t is large. Therefore, we can construct 
integer lattice coverings as nearly eflBcient as desired to a given real lattice covering, thus 
giving asymptotic results for the present case of the degree-diameter problem. The precise 
result is as follows. 

Theorem 9. Let eji be the best possible efficiency for a lattice covering of R*^ by Si, 
and let Szik) be the best possible efficiency for a lattice covering of Z'^ by Sk- Then 
Szik) = £r + 0{k~^). The same applies to S[ and S'^. 

Proof. Let L be a lattice giving a lattice covering of R*^ by 5*1 with efficiency £r. Let 
vi, . . . , v^; and c be as in Theorem 8. Given a large integer k,\eit = k — c, and let L be the 
integer lattice approximating tL as in Theorem 8, generated by wi, . . . , w^. Since L is an 
integer lattice and St+c+L = R'', we have Sk+L = Z^. Let M and M be the d x d matrices 
whose rows are and Wj, respectively; then det M = {2^/d\)sii and detM < \Sk\£z{k)-, 
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which imphes det{k-'^M) < {2'^/dl_+0{k-^))ez{k). But M = tM + 0{l) = kM + 0{l), so 
k-'^M = M+0{k-'^), so det{k-^M) = det M +0{k-^); this imphes ez{k) > eR + 0(A;-^). 

On the other hand, Proposition 4(a) states that any lattice that gives a covering of Z*^ 
by Sk also gives a covering of R'^ by Sk+d/2- If the efficiency of the former is £z(fc), 
then the efficiency of the latter is (|S'fc|/ yoI{S k+d/2))£z{k) = (1 + 0{k~^))ez{k); hence, 
(1 + 0(/c-i))£z(/c) < £R, so ezik)_< £r + 0{k-^). 

The same argument works for S[ and S'f^. ■ 

Again, the same applies to other shapes as well. Combining this with the known sizes 
of the sets Sk and S'^ gives: 

Corollary 10. Let d be a fixed positive integer. 

(a) If e-R, is the best possible efficiency for a lattice covering oflV^ by Si, then the size 
of the largest possible undirected Cayley graph of an Abelian group on d generators with 
diameter at most k is {2'^eji/dl)k'^ + 0(k'^~^). 

(b) If e'^ is the best possible efficiency for a lattice covering oflV^ by S'l, then the size of 
the largest possible directed Cayley graph of an Abelian group on d generators with diameter 
at most k is {eyd\)k'^ + 0{k'^-'^). ■ 

The above assumes that the upper limit £r on the efficiency of a lattice covering of 
by Si is actually attained. To see that this is the case, first note that there is certainly a 
lattice covering with positive efficiency, say Sq. Now consider the possible sets of generating 
vectors vi, . . . , for a lattice L giving a covering of efficiency at least Eq. We may assume 
that, if Vj is one of the generators and w is an integral linear combination of the other 
generators, then \vi\ < |v^ + w|; otherwise, just replace Vj with + w to get a smaller 
set of generating vectors for L, and iterate until no more such reductions are possible. It 
follows that, if, say, vi is the longest of the vectors v^, then the angle between vi and 
the hyperplane P spanned by V2,...,Vd is bounded below by a positive number (tt/S 
for d = 2, somewhat less for higher d). But the distance from vi to P is at most the 
diameter of Si (i.e., 2), since otherwise Si + L would consist of 'hyperplanes' of copies 
of 5^1 with gaps in between. Putting these together, we get a fixed upper bound B on 
the length |vi|, and hence on all of the lengths v^. But we also have a positive lower 
bound b on the determinant det(vi, . . . , v^), namely Sq vo^^i). One can now show that 
there is a fixed number M such that, if mi, . . . , G Z and \mi ! + •••+ \md\ > M, then 
|wivi H — ■ + m^Vdl > 3. Hence, the finitely many translates + mi vi + ■ — (- mdVd with 
|wi I + ■ ■ ■ + |md| < M will have to cover all points at distance < 2 from the origin. Now, the 
set of sequences vi, . . . , with all |vi| < B such that the translates 5'i+miVi + - ■ ■ + mdVrf 
for |mi I + ■ — h \md\ < M cover all points at distance < 2 from is a compact set, so there 
is a sequence vi, . . . , Vd in this set for which det(vi, . . . , Vd) is maximal; this sequence of 
vectors generates a lattice covering of R'' by Si (by Proposition 6) with maximal possible 
efficiency. The same argument works for any compact shape of positive volume, such as S'l. 
(Presumably this argument is well known, but the authors were not able to find a reference 
for it.) 

As we will see in a later section (for the case d = 3, but the argument is general), the 
determination of actual values for en, and for a particular dimension d can in principle 
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be reduced to the solution of a finite number of nonfinear optimization problems (each of 
which requires maximizing a degree-d polynomial function over a region which is a convex 
polytope in R'' ) . Unfortunately, this finite number is extremely large, so the actual values 
are not known for d > 2. (For d = 1 we trivially have e^i = s'^ = 1. For d = 2 we will see 
below that £r, = 1 and e'^ — 2/3.) The computations described in later sections lead us 
to conjecture that, for d = 3, £r = 8/9 and = 63/125. 

Two generators, undirected. We now begin to consider the results that can be obtained 
for specific values of d. As noted previously, the case d = 1 is trivial, so we will start with 
the case of undirected Abelian Cayley graphs on two generators and a given bound k on 
the diameter. The results in this subsection are not new. 

As noted before, the diameter of the ordinary toroidal mesh x is [m./2\ + [n/2j. 
Hence, the largest such mesh with diameter < k is the one where m is k rounded up to the 
nearest odd integer and n is k + 1 rounded up to the nearest odd integer. This corresponds 
to the lattice covering of Z^ by Sk using the lattice L = mZ x nZ; the efficiency of this 
covering is mn/{2k^ + 2k+ 1), which tends to 1/2 as grows large. So one can hope for 
better results. 

To get these, consider the real rotated square Sk- There is obviously a lattice tiling 
using this square; it is just a rotated orthogonal grid with spacing ^/2k between lines. The 
lattice Li for this tiling is generated by the vectors {k,k) and {—k,k); as expected, the 
corresponding determinant is 2A;^, which is equal to the area of Sk- 
it now follows from the approximation results that we can get lattice coverings of Z^ 
using Sk with efficiencies that approach 1 for large k. However, we do not need the general 
results here; since Li is already an integer lattice, we can simply note that Sk + Li = 
implies Sk + Li — Z'^. So this gives a lattice covering using Sk whose index is |Z^ : Li| = 
2A;^, which is better than that from the best toroidal mesh for all A; > 3; for large /c, the 
efficiency approaches 1. 

The corresponding Cayley graph T? jLx turns out to be quite simple to describe. The 
2kxk rectangle {1, . . . , 2k} x {1, . . . , /c} contains exactly one point from each coset of Li, so 
it can serve as a set of vertices for the graph. Adjacent points in the rectangle (horizontally 
and vertically) are connected as in the usual mesh. Horizontally, one has the usual toroidal 
connections at the ends: (1, j) is connected to {2k, j)- But vertically, there is an offset of k: 
(z, 1) is connected to {i + fc, fc) if i < fc, or to (z — /c, fc) if z > k. This is just like a 2k x k 
toroidal mesh, except that the torus is twisted halfway around before the long edges are 
glued together. This twist allows one to double the number of vertices in a k x k toroidal 
mesh while increasing the diameter by at most 1 (there is no increase if k is even). 

A number of the useful properties of ordinary toroidal meshes apply with very little 
change to twisted toroidal meshes. For instance, since the new mesh is still just a rectan- 
gular mesh with extra connections at the boundary, it is easy to map a simple rectangular 
grid into the mesh, by simply ignoring the boundary connections. 

Another nice property of toroidal meshes is that it is easy to find a shortest route from 
one node to another: just check for each coordinate separately which of the two possible 
directions gives a shorter path, and put the results together. Finding optimal routes is 
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only slightly more complicated for the twisted toroidal mesh. To see this, consider the 
given 2k X k rectangle as half of a 2A; x 2k rectangle; each node in the first half has 
a copy {i ± k,j ± k) in the other half. This larger rectangle is then copied periodically 
without further twists to cover Z^; in other words, the 2k x k twisted toroidal mesh is 
just a, 2k X 2k toroidal mesh where and {i ± k,j ± k) are identified as a single node. 
Therefore, to find an optimal route from to in the twisted mesh, apply the 

ordinary 2k x 2k toroidal mesh routing algorithm to find optimal routes from to (i', j') 
and to {i' ± k^j' ± A;), and choose the shorter of the two. 

In the real case the lattice Li gave a perfect tiling of using 5'^, since boundary 
overlap did not count; but in the integer case the boundary overlap reduces the efficiency 
slightly, from 1 to 2/c^/(2/c^ + 2/c + 1). It turns out that if one uses a slightly modified 
lattice, namely the lattice L2 with generating vectors (/c, k + 1) and (— /c — 1, A;), then one 
gets a covering of TP' by copies of with efficiency 1 (i.e., a tiling). See Figure 1. We 
therefore get: 

Theorem 11 (multiple authors). The largest possible size for the undirected Cayley graph 
of an Abelian group on two generators having diameter k is 2k^ + 2A; + 1. ■ 

This result has appeared in various forms in a number of places (usually stated so as to 
apply only to cyclic Cayley graphs, but since the optimal Abelian Cayley graphs turn out 
to be cyclic, the results are basically equivalent). See, for instance, Boesch and Wang [4] 
or Yebra el al. [23]; the Bermond-Comellas-Hsu survey [3] has many additional references. 

The tiling in Figure 1 is in Yebra el al. [23], among other places; it even appears in native 
artwork of the southwestern United States, and may date back to the ancient Aztecs, who 
did use the stepped diamond shape in temple ornamentation. (This shape is now commonly 
known as the Aztec diamond, a term coined by J. Propp.) However, it is unlikely that the 
Aztecs were motivated by the desire to construct efficient parallel computation networks. 

It is easy to see that the lattice tiling of by Sk, or of by the Aztec diamond, 
is unique except for a possible refiection about the line x = y; this just corresponds 
to interchanging the two generators for the Cayley graph. Therefore, the Cayley graph 
attaining the bound in Theorem 11 is unique up to isomorphism. 

Since the point {2k + 1, 1) is in L2, we have e2 + L2 — {—2k — l)(ei + L2) in T? jL^- 
Hence, T? jL^ is a cyclic group, generated by ei + L2 alone. It is isomorphic (not only as 
a group but as a Cayley graph) to Z2fc2_|_2fc+i with the generating set {1,2A;^}. One may 
choose to replace the second generator by its inverse, making the generating set {1, 2A;-|- 1}; 
other generating sets can be used as well. 

For layout purposes, one may just arrange the nodes in the form of the diamond and 
connect the boundary nodes as specified by L2, but it is probably more convenient to use 
the almost- rectangular shape outlined in Figure 1 (a {2k + 1) x rectangle with an extra 
partial row of length A; + 1). The boundary connections are similar to those for the twisted 
toroidal mesh given before, but now there is also a slight twist when connecting the short 
sides; there is a drop of one row when wrapping around from right to left. This layout 
shows that one can embed a rectangular grid into this graph so as to use almost all of the 
nodes. 
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Figure 1. Lattice tiling of T? using S}~ (shown for = 3). 

Two generators, directed. We now describe the largest possible directed Cayley graph 
of an Abelian group on two generators with diameter bounded by /c. As in the preceding 
subsection, the two-generator results here are already known. 

The best toroidal mesh in this case is Zj„ x Zj„/, where m — |_/c/2j +1 and n— [/c/2] +1; 
this gives size mvnl — [(/c + 2)^/4j, which is about 1/2 of j^^l, so one can hope to do better. 

However, one is not going to get perfect efficiency in this case. One can easily tile the 
plane with triangles such as if one is allowed to rotate them, but this is not possible 
using only a lattice of translated copies of a triangle. The exact minimum density for a 
lattice covering of the plane by triangles was computed by I. Fary in 1950; we will give a 
different proof of his result here, and then give the analogue for 7?. 

Theorem 12 (Fary [13]). The minimum density for a lattice covering of by triangles 
is 3/2. Equivalently, the maximum efficiency is 2/3. 

Proof. Since the density and efficiency of a lattice covering are invariant under affine 
transformations, it does not matter which triangle we work with, so, for slight convenience, 
let us work with the isosceles right triangle S[ . 

One can attain the efficiency 2/3 by using the lattice with generating vectors (1/3, 1/3) 
and (2/3, —1/3). This corresponds to a tiling of the plane using an L-tromino that takes 
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Figure 2. A subset of the triangle S[ which tiles the plane. 

up 2/3 of S[, as shown in Figure 2. Or one can cut off all three corners of the triangle to 
get a hexagon that tiles the plane. 

Now suppose that we have a lattice covering of using S[ and the lattice L; we must 
show that the efficiency of the covering is at most 2/3. Let -< be a linear ordering of L 
compatible with addition which is defined by primarily ordering points {x,y) according 
to the sum x + y (so points that are farther out in the direction (1, 1) come later in the 
ordering) and breaking ties (if any) by distance in some other direction. 

Let AB be the hypotenuse of S[. Only finitely many of the L-translates of S[ lie near S[; 
of these, the ones of the form .S^ + v for v must cover the points which are near AB on 
the side away from S[. Since the union of finitely many translates of S[ is closed, AB itself 
is covered by finitely many translates + v with v :^ 0. Find such a covering of AB with 
as few translates as possible, say S[ + Vi, . . . , S'l + v^, where Vj >- 0. 

Note that, since S[ + Vi must intersect AB, it contains one of the endpoints A and B 
if and only if the coordinates of are not both positive. We may assume that at most 
one of the vectors has both coordinates positive. For if there are two such, let them 
be S[ + Vi and S[ + vj where -< Vj. Then — :^ 0; since Vj has both coordinates 
positive, we have 

{S[ + Vj) nABC(S[+ Vj - Vi) n AB. 

Furthermore, Vj — Vj cannot have both coordinates positive, because, if it did, we would 
have 

(S[ + Vj) n AB C (^( + Vi) n AB, 

so S[+Vj would not have been needed in the covering of AB, contradicting the minimality 
of m. So we can replace S[ + Vj with S'l + Vj — Vi to get another covering of AB using 
fewer vectors with both coordinates positive. Repeat this until only one such vector is left. 

Since each translate S'l + Vj is convex, its intersection with AB is a segment or a 
point. Therefore, at most one of these translates can contain A, since otherwise one of the 
intersections + Vj) fl AB would include another such intersection, making the latter 
translate superfluous and contradicting the minimality of m. Similarly, at most one of the 
translates S[ + Vj contains B. Putting these facts together, we conclude that we needed 
at most three of the translates -S^ + v with v :^ to cover the segment AB. 
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In other words, there are points P and Q on AB such that each of the three segments 
AP, PQ, and QB is covered by one of the translates S'i + v with v :^ 0. Let h,l2, h be 
the lengths of these three segments; then + ^2 + ^3 = V^- Note that, if S'J + v covers AP, 
then 5'( + V covers the entire isosceles right triangle below AP whose hypotenuse is AP; 
the area of this triangle is if/ 4:. Similar statements hold for PQ and QB. So we have three 
disjoint triangles included in S[ which are covered by translates + v with v 0, and 
the total area of these triangles is {if + I2 + ^3)/4. 

By the proof of (the real version of) Lemma 3, if we let T be the part of §[ which 
is not covered by any translate + v with v 0, then T gives a lattice tiling of 
using L, so the efficiency of the covering using S[ and L is Area(T)/Area(S'(). We have 
Area(S'i) = 1/2 and Area(r) < 1/2- {ll + ll+lD/A. A standard minimization shows that, 
if 1-^ + 12 + 13 = V2, then + ^2 + > 2/3 (with equality only when li = I2 = Is = V2/3). 
Therefore, Area(T) < 1/3, so the efficiency of the covering by S[ and L is at most 2/3. ■ 

It now follows from Theorem 9 that the largest possible index |Z^ : L| for an integer 
lattice L giving a lattice covering of by S'j^ is approximately (2/3)|5'^|, or about A;^/3, 
for large k. However, we can actually get an exact answer rather than an approximation. 

Theorem 13 (mainly Wong and Coppersmith [22]). The largest possible index \T? : L\ 
for a lattice L giving a lattice covering of T? by S'/^ is [{k + 2)^/3j . 

Proof. We will give a discrete form of the proof of Theorem 12. Let a be {k + 2)/3 rounded 
to the nearest integer, and let 6 = A; + 2 — 2a (so 6 is also about (fc + 2)/3). Let Tk be the 
set of (i,i) in such that z, j > 0, min(i, j) < a, and max(i, j) < a + b. Then Tk C S'j^, 
since any {i,j) in T^ satisfies i + j <a — 1 + a + b — 1 = k. The set T^ looks like the 
L-tromino from Figure 2, and it tiles Z^ using the lattice with generating vectors (a, a) 
and (a + b, —b). Therefore, this lattice gives a covering of Z^ using 5'^, and its index is 
a(a + 2b), which works out to be [{k + 2)^/3] . 

Now, suppose we have a lattice covering of Z^ using S'j^ and a lattice L; we must show 
that |Z^ : L\ < [{k + 2)^/3j. Define the linear order -< of L as before. Let A and B be 
the points (0, k + 1) and {k + 1, 0); then the segment AB contains k + 2 integral points, 
which must be covered by translates S'^ + v where v e L and v :^ 0. 

Let vi , . . . , be a list of as few vectors as possible in L such that Vj :^ and the 
translates •S'^^^^+v^ of S'j^_^-^ cover all of the integral points on AB. Then the same argument 
as for Theorem 12 shows that m is at most 3. Hence, AB can be broken up into three 
segments AP, P'Q, and Q'B (where P and P' are adjacent integral points on AB, as are 
Q and Q'), each of whose integral points are covered by one of the translates 'S'^+i + v 
with V >- 0. Let h,l2,h be the numbers of integral points on these segments; then 
h + h + h = k + 2. 

If 5*^+1 + V covers the integral points on AP, then it covers all of the integral points in 
the isosceles right triangle below AP and having AP as its hypotenuse. In fact, all of these 
points other than those on AP itself are covered by 5^ + v; there are (/f — /i ) /2 such points, 
and they are all in S'f.. Similarly, the segments P'Q and Q'B give (/§ — ^2)/2 + (if — l3)/2 
more points of S'f. which are covered by translates + v with v :^ 0. 
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As in the proof of Lemma 3, let T be the set of points in S'^ that are not in 5^ + v for 
any v :^ 0; then |T| = : L|. The calculations above show that 

|T| < l-S^I + (li + ^2 + h)/2 - (ll + /2 + ll)/2. 

Here \S'j^ \ — {k + l){k + 2)/2, and I1+I2 + h is just k + 2. Given h + I2 + h, we minimize 
II + I2 + ^3 by making the numbers Zi, Z2, ^3 as close to equal as possible; in this case, this 
means that the minimum occurs when two of them are a and the third is b. Therefore, 

{k + l){k + 2) k + 2 _ 2a^ + 6^ 
1^1- 2 +2 2 ' 

which simplifies to \T\ < [{k + 2)^/3j, as desired. ■ 

Corollary 14 (mainly Wong and Coppersmith [22] ) . The largest possible size for the di- 
rected Cayley graph of an Abelian group on two generators having diameter k is [(A;+2)^/3j . 

■ 

Again it is hard to be historically accurate here, because different authors have pre- 
sented results in quite different ways; see the Bermond-Comellas-Hsu survey [3] for more 
information and references. 

Let a, 6, and be as in the proof of Theorem 13; then the set gives a suitable 
layout for a network realizing this Cayley graph. In addition to the mesh connections from 
(i, j) to {i + and + 1) within T^, one will also need wraparound connections from 
(a + 6 — 1, j) to (0, j + b) for j < a, from (a — 1, j + a) to (0, j) for j < b, from {i,a + b — 1) 
to {i + b, 0) for i < a, and from (i + a, a — 1) to (i, 0) for i < b. 

In the case a = b, one can use an alternate layout in the form of a 3a x a rectangle, 
with wraparound connections from {3a, j) to (1, j) and from {i,a) to {{i + a) mod a, 1). 
This is just a variant of the twisted toroidal mesh where the long dimension is twisted 
by a factor of 1/3 rather than 1/2; it is convenient for construction and for embedding 
a rectangular grid without boundary connections into the network (although this is not 
particularly useful in the directed case). li a ^ b, then one gets a rectangle with some 
missing nodes or extra nodes along part of one edge, and the cross-connections are slightly 
messier. 

If k = 1 (mod 3), so that a = b = {k + 2)/3, then one can see from the proof of 
Theorem 13 that the lattice L with generating vectors (a, a) and (a + 6, —6) is the unique 
lattice attaining the bound in the theorem, and hence the Cayley graph attaining the bound 
in Corollary 14 is also unique. However, if A; ^ 1 (mod 3), so that a and b diff'er by 1, 
then there arc two more lattices attaining the bound: the lattice L with generating vectors 
(a, b) and (2a, —a), and the mirror image with generating vectors (6, a) and (—a, 2a). The 
latter two give Cayley graphs that are isomorphic to each other, but not to the Cayley 
graph of Z^/L (if A; > 1), because the Cayley graph of Z^/L has cycles of length 2a while 
that of Z^/L does not. Therefore, if A; > and k^ 1 (mod 3), then there are exactly two 
Cayley graphs meeting the bound of Corollary 14. 
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If A; > 1 and k = 1 (mod 3), then the optimal group Z^/L is not cychc; it is isomorphic 
to Z30 X Za by an isomorphism sending ei and 62 to (1,0) and (3a — 1,1). On the other 

hand, if k ^ 1 (mod 3), then (2a + &, a — b) is in L and a — b = ±1, so 62 is a multiple of ei 
in I? jL and hence T? jL is cyclic. Similarly, I? jL is cyclic, since (3a, b — a) G L. One can 
get the corresponding Cayley graphs directly from the cyclic group Z\_(^k^2y /3} by using 
the generator pairs {1, (2a + b)/{b — a)} and {1, 3a/ (a — 6)}, respectively. 

Three generators, undirected. For d = 3, we must consider three-dimensional lattice 
tilings by the regular octahedron and its discrete approximation S^- These shapes do 
not tile space perfectly, and the best possible efficiency for a lattice covering of by 
appears to be still open (although there is a good guess, as we shall see). So we will apply 
our results in reverse, using computed results about the degree-diameter problem to obtain 
information about lattice tilings by octahedra. 

The best three-dimensional toroidal mesh with diameter k is Z260+1 x Z261+1 x Z262-1-1, 
where hi = [{k + i)/3j; this has about (8/27)/c^ vertices for large k. This corresponds to 
the covering of by Si using the cubic lattice (2/3)Z^; this covering has efficiency 2/9. 

It turns out that a good lattice to use for coverings with regular octahedra is the body- 
centered cubic lattice, defined most simply as the set L^cc of points in Z^ whose coordinates 
are all odd or all even. If x is an arbitrary point of R^, then x lies in or on one of the unit 
cubes with vertices in Z^; two opposite corners of this cube will be in Lbcc, say v and w. 
Then each coordinate of x lies between (inclusively) the corresponding coordinates of v 
and w, so, letting S be the metric on R^, we have 5(v, x)+5(x, w) = (5(v, w) = 3. Hence, 
either 5(v, x) < 3/2 or 5(w,x) < 3/2. This shows that |S^3/2| + -^bcc = R*^- Now, Lbcc has 
generators (2,0,0), (0,2,0), and (1,1,1), giving a matrix with determinant 4, while the 
volume of 53/2 is 9/2, so this lattice covering of R"^ has efficiency 8/9. A fundamental 
region for the lattice can be obtained by truncating each of the corners of the octahedron, 
giving an Archimedean solid whose faces are eight regular hexagons and six squares. 

The same reasoning shows that, for > 1, one can get a lattice covering of Z^ by 5"^ 
using the slightly distorted body-centered cubic lattice Lbcc(fti, ^2, ^3) with generating 
vectors (2ai,0,0), (0,2a2,0), and (ai,a2,a3), where a^ = [{2k + This gives a Cayley 

graph of size 4aia2a3, or approximately (32/27)fc^ for large k. This is an improvement 
over the best toroidal mesh of diameter k; it is about 4 times as good for large k. 

One can lay out the Cayley graph for Z^/Lbcc(oi, a2, a3) in the form of a 2ai x 2a2 x a3 
mesh. Opposite 2ai x as sides are connected to each other as in the usual toroidal mesh, 
but the toroidal connections between the top and bottom 2ai x 2a2 sides are twisted in two 
directions: node (ii,i2,«3) is connected to node (ji ± ai, j2 ± «2, 1), where the signs are 
chosen to give numbers between 1 and 2aj, inclusive. Routing algorithms and embeddings 
of rectangular grids work here just as they did in the two-dimensional version. 

Two questions now arise. First, can one improve the efficiency by making small adjust- 
ments to the discrete lattice, as we did in the two-generator cases? Second, can one get 
better results by using a completely different lattice? The answers to these questions are 
not immediately clear, so we will approach the problem from another direction. 

One can write computer programs to examine various groups, choose all (or at least 
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many) possible sets of a certain number of generators for the group, and compute the 
resulting diameters. M. Dinneen has performed many such computations, some using 
exhaustive search of generator sets and others using random sampling, on a number of 
different kinds of groups, resulting in new best-known graphs for the degree-diameter 
problem; see, for instance, Dinneen and Hafner [10]. Some of Dinneen's earlier unpublished 
computations were for Abelian (usually cyclic) groups of diameter up to 10 on various 
numbers of generators. 

The authors have written a program to extend these calculations. The program does 
an exhaustive search of generating sets for each Abelian group, but avoids examining 
many generating sets which give Cayley graphs isomorphic to ones already examined; 
for instance, in the case of a cyclic group Z„, one may assume that the first generator 
is a divisor of n. Here 'exhaustive search' means that all Abelian groups of size up to 
jS'fcl = {4k^ +6k'^ +8k+3) /3 were examined, so the results definitely give the largest possible 
Cayley graph of an Abelian group with diameter k. The program uses bit manipulations 
adapted from (but simpler than) those of Dougherty and Janwa [11], which gave algorithms 
for diameter computations for Cayley graphs of Abelian groups of exponent 2. 

It turns out that, for each k for which the calculation has been done so far (up to 14), 
the best Abelian Cayley graph has been obtained from a cyclic group. The results of the 
computation are shown in Table 1. This extends (and corrects an erroneous final entry in) 
a similar table given by Chen and Jia [5] . 
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25 


9 
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6 


377 
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320 


333 


1,9,73 


.883289 


.592000 


7 
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175 


500 


515 


1,46,56 


.895652 


.628944 


8 


833 


245 


720 


737 


1,11,133 


.884754 


.644700 


9 


1159 


343 


1008 


1027 


1,13,157 


.886109 


.665371 


10 


1561 


441 


1372 


1393 


1,92,106 


.892377 


.686940 


11 


2047 


567 


1792 


1815 


1,15,241 


.886663 


.696960 


12 


2625 


729 


2304 


2329 


1,17,273 


.887238 


.709953 


13 


3303 


891 


2916 


2943 


1,154,172 


.891008 


.724015 


14 


4089 


1089 


3600 


3629 


1,19,381 


.887503 


.730892 



Table 1. Best undirected Cayley graphs of cyclic groups, three generators. 

The first column is the desired diameter k. The second column gives the largest size one 
could hope for of an undirected Cayley graph of an Abelian group on 3 generators. The 
next two columns give the sizes attained by the best possible ordinary toroidal mesh and 
the twisted toroidal mesh described above. Next comes tt-c, the computed largest n such 
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that Tin has three generators giving it an undirected diameter of k. Then comes a triple of 
generators of Z„_, attaining this diameter (this is not always unique, but only one generator 
set is given here). The final two columns gives the efficiencies of the corresponding lattice 
coverings of T? by Sk and of R"^ by >Sa;+3/2 (see Proposition 4). 

Some interesting observations can be made from Table 1. First, note that the twisted 
toroidal meshes do almost as well as the optimal cyclic groups. Also note that the numbers 
in the second-to-last column do seem to be getting close to 8/9 for larger A;; this provides 
evidence that the body-centered cubic lattice gives the best lattice covering of by S\ . 

One can confirm this more strongly by reconstructing the lattices L for which T? jL gives 
these optimal cyclic groups. For instance, look bX k — 10, for which we have the cyclic 
group Zi393 with generating set {1,92, 106}. There is a unique homomorphism from 
to Zi393 which sends ei, 62, 63 to 1, 92, 106, and the desired lattice L is just the kernel of 
this homomorphism; this means that 

L = {(xi,X2,X3) e Z^: xi + 92x2 + 106x3 = (mod 1393)}. 

One can easily find three vectors in L, namely (1393,0,0), (92,-1,0), and (106,0,-1); 
the matrix with these three vectors as rows has determinant 1393 = : L|, so these 
vectors generate L. Now one can perform elementary operations to reduce these vectors to 
a smaller set of generators for L, such as (7, 7, 7), (8, —7, 6), and (6, 8, —7). These vectors 
are quite close to the vectors (7, 7, 7), (7, —7, 7), and (7, 7, —7), which generate a scaled-up 
body-centered cubic lattice (in fact, the latter lattice gives the twisted toroidal mesh of 
size 1372 mentioned in the table). Similarly, one finds that the other lattices corresponding 
to the generators in Table 1 are almost body-centered cubic. 

There are definite patterns in Table 1; every third k gives groups and generators of the 
same form. These patterns can be generalized, giving the following result. 

Theorem 15. For all k > 0, there is an undirected Cayley graph on three generators of 
an Abelian (in fact, cyclic) group which has diameter k and size n, where 

{(32A;3 + 48A;2 + 54/c + 27) /27 z/ /c = (mod 3), 
(32A;3 + 48A;2 + 78A; + 31) /27 ifk=l (mod 3), 
(32A;3 + 48A;2 + 54A; + ll)/27 ifk = 2 (mod 3). 



Proof. We will show the existence of lattices C 7? such that /Lj^ is cyclic, Sk + Lk = 
Z'^, and IZ"^ : L| is the n specified in the theorem. 

Let a = [2A;/3]. For each k, we define by specifying three generating vectors 
vi, V2, V3 for it, as follows: 

{{a+l,a,a),{a,—a,a+l),{a+l,a—l,—a—l) if A; = (mod 3), 

{a, a, a), {a+1,— a, a— 1), {a— 1, a+1,— a) if A; = 1 (mod 3), 

{a, a, a— 1), {a— 1,— a, a), {a, a— 1,— a) if A; = 2 (mod 3). 
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A simple determinant computation shows that |Z^ : ife| is (2a^ + a+l)(2a+l), Aa^ + Sa, 
or (2a^ — a + l)(2a — 1) in the respective cases A; = 0, A; = 1, or A; = 2 (mod 3). Since 
a is respectively 2A;/3, (2A; + l)/3, or (2A; + 2)/3, the index |Z^ : Lk\ works out to be the 
desired value n. 

For A; = (mod 3), the following vectors are in L^: 

V2 + V3 = (2a+l,-l,0), 
vi + (2a-l)v2 + 2av3 = (4a^+2a+l, 0, -1). 

Hence, we have e2 = (2a+l)ei and 63 = (4a^+2a+l)ei in Z^/Lk, so ei generates 1? jL^. 
So 7? / Lk is isomorphic to Z^, via an isomorphism taking ei, 62, 63 to 1, 2a+l, 4a^+2a+l. 
Similarly, for A; = 1 (mod 3) we have 

av2 + (a-l)v3 = (2a^-a+l, -1, 0), 
(a+l)v2 + av3 = (2a2+a+l, 0, -1), 

so Z^/Lfc is isomorphic to Z„ with generators 1, 2a^— a+1, 2a^-|-a+l; and for A; = 2 
(mod 3) we have 

V2 + V3 = (2a-l,-l,0), 
vi + (2a-l)v2 + 2av3 = (4a^-2a+l,0, -1), 

so Z^/Lfc is isomorphic to Z^ with generators 1, 2a— 1, 4a^— 2a+l. 

It remains to show that Sk + Lk = ■ We will do only the case k = 1 (mod 3) here; the 
other two cases are handled by the same method, but with a few more subcases because 
of less symmetry. 

For A; = 1 one just has to show that Z7 with generators 1, 2, 4 has diameter 1, and this 
is trivial to do directly; so we may assume A; > 1 and hence a > 1. 

Let V4 = vi — V2 — V3 = (—a, a— 1, a+1). Then the vectors iv^ for z = 1, 2, 3, 4 give one 
member of strictly within each of the eight octants of Z^, and all of the coordinates of 
these vectors have absolute value at most a+1. 

We must show that each x e Z^ is in Sk + Lk- This is equivalent to showing that there 
is a member w of Lk such that x — w G 5"^, which in turn is equivalent to 5(x, w) < k, 
where d is the (Manhattan) metric on Z^. Note that if x, y, z are such that each 
coordinate of y is between (inclusively) the corresponding coordinates of x and z, then 
5(x, y) + 5(y, z) = 5(x, z). Prom now on, we will state this situation more briefly as "y lies 
between x and z." 

Suppose we are given x G Z^. The idea is to repeatedly reduce x by adding members 
of Lk to it, until one reaches a vector which is within /^-distance A; of or some other 
known member of Lk- 

The first thing we will do is reduce x to a vector whose coordinates all have absolute 
value at most a+1. Suppose x does not already have this property. Let v be one of 
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the vectors ±Vj {i = 1,2,3,4) such that the coordinates of v have the same signs as the 
corresponding coordinates of x; if a coordinate of x is 0, then either sign is allowed for the 
corresponding coordinate of v. Now look at x' = x — v. If a coordinate of x has absolute 
value < a+1, then the corresponding coordinate of x' will also have absolute value < a+1, 
because of the sign matching and the fact that the coordinates of v have absolute value 
< a+1. If a coordinate of x has absolute value > a+1, then the corresponding coordinate 
of x' will be strictly smaller in absolute value. Therefore, repeating this procedure will 
lead after finitely many steps to a vector whose coordinates all have absolute value at most 
a+1. 

If this new x lies between and one of the vectors iv^, then we have 5(0, x)+5(x, ±Vi) = 
(5(0, ±Vi). But all of the vectors ±Vi satisfy 5(0, iv^) — 2k + l; since 5(0, x) and 5(x, iv^) 
are both integers, one of them must be at most k, so we are done with this x. 

We now break into cases depending on which octant the new x lies in. Since Lk is 
centrosymmetric, we only need to handle the octants containing vi, V2, V3, and V4. Also, 
Lk is invariant under cyclic permutations of the three coordinates, since these leave vi 
fixed and permute V2, V3, V4; hence, we may assume that the new x is in the octant of vi 
or the octant of V2. 

First suppose that x is now in the octant of Vi (all three coordinates nonnegative). If 
X is between and vi, we are done. If two or more of the coordinates of x are equal to 
a+1, say (by cyclic symmetry of Lk) x = (a+1, a+1, r), then we have 5(x, vi) < k unless 
k = A and r = 0, in which case 5(x, vi + V3) = 4 = k. 

If X has exactly one coordinate equal to a+1, say x = (a+1, r, s) with < r, s < a, then 
we can subtract vi from x to get x' = (l,r— a, a), which is in the octant containing 
— V4. If x' lies between and — V4, we are done. If not, then r = 0. Now let x" = x' + V4 = 
(—a+1, —1, s+1), which lies between and — V3 unless s = a, in which case x = (a+1, 0, a) 
and 5(x, V2) = a+1 < k. 

The procedure when x is in the octant of V2 is similar. Either x = (r, s, t) lies between 
and V2, or s — — a— 1, or t > a. In the latter cases, let x' = x — V2. When s = —a—1, we 
have x' = (r— a— 1, — l,t— a+1); either this lies between and one of the vectors ±v^, or 
X + V3 does. When s > —a but t > a, try x' — V4; either it lies between and some ±Vj, 
or it is (—1, 1, —a), (a, 1, —a), or (a, 1, —a+1). These last three lie within 5-distance k of 
V3 — Vi, V3, and either V3 + V2 — Vi or V3 + V2, respectively. ■ 

The authors conjecture that the graphs given by this theorem are actually the largest 
undirected Cayley graphs of Abelian groups on three generators for each diameter k. 

This conjecture would imply that the lattice covering of R"^ by 53/2 using the lattice L^cc 
is optimal; that is, 8/9 is the best possible efficiency for a lattice covering by regular 
octahedra. The latter statement seems quite plausible, but remains unproved at this 
point. We can prove a partial result, though, that a "small" adjustment to Lbcc cannot 
improve the covering: 

Theorem 16. Among those lattices L for which + L = R^, the lattice Lbcc is locally 
optimal; that is, for any other lattice L sufficiently near Lbcc such that S3/2 + L = R^, the 
efficiency of the covering using L is less than 8/9. 
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Proof. Let us use the vectors vi = (1, 1, 1), V2 = (1, —1, 1), and V3 = (1, 1, —1) as gen- 
erating vectors for Lbcci then a nearby lattice L will be generated by nearby vectors 
= (ai,&i,ci), w'2 — (02,62,02), and Vg = (a3,^>3,C3). Wc can concatenate the three 
vectors v']^,V2,V3 to get a single vector v' in R^; similarly, let v be the concatenation 
vi,V2,V3. Let F{y') be the determinant of the matrix with rows vi,V2,V3. Note that 
F(v) = 4; we must see that this point is a strict local maximum of -F(v') for those points v' 
satisfying the constraint that §^12 + L = R^. We compute that the gradient of F at the 
point V is g = (0, 2, 2, 2, -2, 0, 2, 0, -2). 

Using the lattice Lbcc, the point (1/2, 1/2, 1/2), in the center of a face of 5^3/2, is covered 
by only two copies of namely S3/2 itself and S3/2 + vi, and it is on the boundary (a 
face) of each of these copies. If the lattice is altered slightly so that these two copies no 
longer touch, then the points in between will not be covered by any copy. In particular, if 
L is near Lbcc but ai + bi + ci > 3, then the point (1/2, 1/2, 1/2 + e) for small positive s 
will not be in ^'3/2 + L. So the constraint 5*3/2 + = R^ gives us the linear inequality 
ai + bi + ci < 3. We will rewrite this as 

ui • v' < 3, where Ui = (1, 1, 1, 0, 0, 0, 0, 0, 0). 

The same argument for points on the other faces of the octahedron gives inequalities 

U2 ■ v' < 3, where U2 = (0, 0, 0, 1, -1, 1, 0, 0, 0), 
U3 • v' < 3, where U3 = (0, 0, 0, 0, 0, 0, 1, 1, -1), 
U4 • v' < 3, where U4 = (—1, 1, 1, 1, —1, —1, 1, —1, —1). 

Next, consider the point (1,0, 1/2). This is in ^3/2 + y for four members y of Lbcc, 
namely 0, vi, V2, and V2 + V3, and it is an edge point of each of these four copies. If the 
lattice is altered slightly, then a gap can open up near this point even if there are no gaps 
between octahedra adjacent at a face as above. 

Specifically, if v' is near v, £ is a very small positive number, and we define the point x 
by the linear equations 

X • (1, -1, -l) = v; -(1,-1,-1)+ 3/2 + £, 

X • (-1, -1, 1) = (v^ + v^) ■ (-1, -1,1)+ 3/2 + e, and 

X • (-1, 1, -1) = v'2 ■ (-1, 1, -1) + 3/2 + e, 

then X will be a point near (1, 0, 1/2) which is not in S3/2 + y for y G {v^, V2, V2 + V3}. 
Adding up the three given equations yields 

X • (-1, -1, -1) = vi • (1, -1, -1) + • (-2, 0, 0) + • (-1, -1, 1) + 9/2 + 3e. 

If the right hand side of this equation is less than —3/2, then x will not be in ^3/2 either, 
and hence will not be in + L. Since s can be arbitrarily small, in order to have 
53/2 + L = R^, it is necessary to have 

v'l • (1, -1, -1) + V2 ■ (-2, 0, 0) + • (-1, -1, 1) > -6. 
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This can be rewritten as 

Us • v' < 6, where U5 = (-1, 1, 1, 2, 0, 0, 1, 1, -1). 

The same argument can be performed using the octahcdra around (1, 0, 1/2) in the opposite 
order, and there are 23 other points on the edges of 6*3/2 where the same configuration 
occurs. But one only gets six distinct inequahties from this; the other five are: 

U6 • v' < 6, where ue = (0, 2, 0, 1, —1, 1, 1, —1, —1), 

U7 • v' < 6, where U7 = (-1, 1, 1, 1, -1, 1, 2, 0, 0), 

Us • v' < 6, where Ug = (1, 1, 1, 0, -2, 0, 1, -1, -1), 

Ug • v' < 6, where Ug = (0, 0, 2, 1, —1, —1, 1, 1, —1), 

Uio • v' < 6, where Uio = (1, 1, 1, 1, —1, —1, 0, 0, —2). 

Note that all ten of these inequahties are satisfied with equahty when v' = v. Hence, 
they can be rewritten as Uj • (v' — v) < for z = 1, 2, . . . , 10. 

One can easily check that the vectors Ui, . . . , U7 are linearly independent; their common 
null space (i.e., the set of w such that Uj • w = for all z < 7) is generated by the 
independent vectors Wi = (1, 0, —1, 1, 0, —1, 1, 0, 1) and W2 = (—1, 1, 0, —1, —1, 0, —1, 1, 0). 
Also, we have 



g = Ui + U2 + U3 + U4 = U5 + Ug = U6 + Ug = Uy + Uiq. 



Let C be the closed cone consisting of all vectors t in the subspace spanned by ui , . . . , U7 
such that Uj • t < for all i < 10. Then the above equations imply that g ■ t < for all 
t in C, and equality can hold only when t = 0. In particular, we have g • to < for any 
unit vector to in C. The set of such to is closed and bounded, hence compact, so there 
is a positive number e such that g • to < —s for all such to. It follows that there is a 
neighborhood t/ of g such that, for any g' in U and any unit vector to in C, g' • to < 0. 
Since C is a cone, we have g' • t < for all g' G and all nonzero t e C. 

We can compute that, for any real numbers r and s, the deteriminant for the lattice 
given by V + rwi + SW2 is 

F(v + rwi + SW2) = 4(1 - r)(l + s)(l + r - s). 

If |r| + |s| < 1, then the numbers 1— r, 1+s, and 1+r—s are positive numbers with arithmetic 
mean 1, so their geometric mean is at most 1; this means that F(v + rwi + SW2) < 4. 
Equality holds only when the above three numbers are equal, which is when r = s = 0. 

Let U' be a convex neighborhood of v so small that (gradF)(v') G U for all v' G U' . 
Now, any vector v' sufficiently close to v can be expressed as v + ti + t2 where ti is a 
(small) linear combination of Wi and W2, t2 is a linear combination of Ui, . . . , Uy, and both 
V + ti and V + ti + t2 are in U'. If v' satisfies the condition 8^/2 + L = and is near 
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to V, then we must have • (v' — v) < for all i < 10, so ■ t2 < for all i < 10 (since 
• ti = 0), so t2 e C. We have F(v + ti) < 4, with equality holding only when ti = 0. 
If t2 is nonzero, then for any t on the segment from v + ti to v + ti + 12 we have t G U' , 
so (gradF)(t) e U, so (gradF)(t) ■ ts < 0; it follows that F(v + ti + t2) < F(v + ti). 
Therefore, -F(v') < F{v), with equality holding only when v' = v. So v gives a local 
maximum of F, as desired. ■ 

It is still possible (though very unlikely) that a lattice quite different from L^cc gives a 
more efficient covering. Theoretically, the search for an optimal lattice can be set up as a 
large optimization problem and solved once and for all, but this appears to be a formidable 
task. 

One could begin this task by considering an arbitrary lattice L such that 3^/2 + L = 
and this covering is reasonably efficient (at least as efficient as the covering from Lbcc)- Such 
a lattice is generated by vectors vi,V2,V3, and we can carefully choose these generators 
so as to limit their lengths. In particular, we can choose Vi to be a nonzero member of L 
with minimal length. We can then choose V2 in L whose distance from the subspace of 
spanned by vi is as small as possible (but nonzero), and adjust V2 by subtracting off an 
integer multiple of vi so as to ensure that the closest integer multiple of vi to V2 is 0. One 
can similarly choose V3 to be as close as possible to (but not in) the subspace spanned by 
Vi and V2. These three chosen vectors will be a set of generating vectors for L. In order 
to have ^3/2 + ^ = R^, it is necessary that the length of Vi be no more than the diameter 
of 5'3/2; there are similar but slightly larger bounds on the lengths of V2 and V3. This 
limits our search for vi,V2,V3 to a compact subset of nine-dimensional space. We must 
find the point in this subset which maximizes det(vi, V2, V3) subject to the constraint that 
^3/2 + L = R3. 

This constraint looks infinitary, but it can actually be reduced to finitely many sets of 
linear inequalities. To see this, note that, using the above upper bounds on the lengths 
of the vectors along with the assumed lower bound on the lattice determinant (the 
covering must be at least as efficient as that from Lbcc); we can get lower bounds on the 
lengths of the Vj, the angles between them, and associated quantities such as the distance 
from V3 to the plane spanned by Vi and V2 . These will allow us to get upper bounds on the 
absolute values of integers ai, 02, 03 such that <53/2 + aiVi -|-a2V2 + a3V3 overlaps or almost 
touches 5*3/2 • (In other words, we get an upper bound on the number M from the proof of 
Proposition 7.) So we only have to consider finitely many of the lattice translates of 5^3/2 
when trying to cover the space near Ss/2 (which is all that is needed, by Proposition 6). 

There are only finitely many configurations (specifications of arrangements and overlaps) 
for these finitely many translates of 5*3/2- For each such configuration, the assertion that 
there are no 'gaps' in the coverage of the space near 5^3/2 becomes a list of linear inequalities 
like the inequalities ■ v' < 6 from the proof of Theorem 16. So we need to optimize a 
cubic function (the lattice determinant) subject to a list of linear inequalities in order to 
find the optimal version of each configuration, and then compare the resulting values to 
find the best configuration. 

Unfortunately, there is a very large number of possible configurations (for an example 
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of the possibilities for complicated configurations, see Figure 4 later in this paper), so this 
finite computation appears to be beyond reach at present. Of course, a different approach 
to the problem might lead to a more feasible computation. 

One might hope to be able to use the arguments of Proposition 7 and Theorem 8 
in reverse, to get an upper bound on the efficiency of lattice coverings of R3 by the 
octahedron Si hy showing that any extremely efficient real lattice covering would lead to 
integer lattice coverings more efficient than what the computation actually found. To do 
this, one would fix a value for the distance p from Proposition 7, and then use the method 
described above to get an upper bound on the number M from that Proposition. If there 
is actually a lattice covering of R"^ by Si using the lattice L generated by vi, V2, V3 having 
a specified large determinant (equivalently, a specified large efficiency), then we can roinid 
the coordinates of these vectors to the nearest multiples of 1/k to get vectors , V2 , Vg 
generating a lattice L'. By Proposition 7, if l/(2fc) is less than 77 = p/M, then we will have 
Si+3p+L' = R'^, so Sk+3kp + kL' = R^. But kL' is an integer lattice; if n is its determinant, 
then this lattice covering will yield an Abelian Cayley graph on three generators with size n 
and diameter at most k + 3kp. The fact that L' is close to L means that we can get a 
lower bound on n from the determinant of L. If the actual computational search showed 
that there is no Abelian Cayley graph of such a size for this diameter, then our original 
assumption that there was a lattice L giving a covering of that efficiency must have been 
false. 

Unfortunately, the constants involved are such that even the large computation done 
so far does not suffice to get a bound less than 1 for the efficiency of L (even if we are 
optimistic enough to assume that M is as small as 3 or 4) . It probably requires searches for 
values of k larger than 500 in order to get actual results from this method; such searches 
are completely out of range at the moment. 

Three generators, directed. For the directed case of three generators, we want to study 
lattice coverings of R^ by the trirect angular tetrahedron S^.. (Since lattice covering effi- 
ciency is affine invariant, it makes no difference which particular tetrahedron we consider.) 
One hopes that one can discretize these coverings to get good lattice coverings of by S"^, 
and hence good directed Cayley graphs. 

The best three-dimensional directed toroidal mesh with diameter k is Zb^+i x Zfo^+i x 
Zfoj+i, where bi — \ i.k+i)/3\; this has about (l/27)/c^ vertices for large k. This corresponds 
to the covering of R"^ by S'l using the cubic lattice (1/3)Z^; this covering has efficiency 2/9. 

It is more difficult to find a candidate for a good covering lattice (or, equivalently, a 
large subset which gives a lattice tiling) for the tetrahedron than it was for the octahedron. 
One possible method is to try to find the three-dimensional analogue of the L-tromino used 
for the triangle; this leads one to consider the tctracube shown on the left of Figure 3. In 
order for the shape to fit into the cdgc-lcngth of the subcubes should be 1/4. It is 
easy to sec that this shape docs indeed tile space, using the lattice generated by (1/2, 0, 0), 
(0, 1/2, 0), and (1/4, 1/4, -1/4) (this is just (l/4)Lbcc); since the shape has volume 1/16 
while S'l has volume 1/6, we get a lattice covering of R^ by S'l with efficiency 3/8. 

The discrete form of this shape, scaled by a factor Si in the z'th dimension, is a subset 
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Figure 3. Two subsets of the tetrahedron S[ which tile space. 

of of size 4S1S2S3 which gives a lattice tiling of Z"^; this subset is included in S'j^, where 
A; = si + S2 + S3 + max(si, .S2, .S3) — 3. Optimizing this for a given k > 1 gives a subset 
of S'l^ which tiles and has size Aa^a^a^, where = + i)/4j. 

One can obtain another lattice covering of by tetrahedra as follows. If one cuts 
off the four corners of a regular tetrahedron at planes passing through the midpoints of 
the edges (so one removes four half-size regular tetrahedra), then what is left is a regular 
octahedron with volume 1/2 that of the tetrahedron. We have a lattice giving a covering 
of by this octahedron with efficiency 8/9; the same lattice therefore gives a covering 
of R^ by the original tetrahedron with efficiency 4/9. 

If one uses an affine transformation to change the regular tetrahedron to the tetrahedron 
S[, then the corresponding lattice will be generated by (1/6, 1/6, 1/6), (1/6, —1/2, 1/6), and 
(1/6, 1/6, —1/2). One fundamental region for this lattice is an affinely distorted truncated 
octahedron. Another one can be obtained by the method of Lemma 3, using an ordering -< 
which orders vectors primarily by the sum of their coordinates; the resulting region is shown 
at the right of Figure 3. This shape consists of 16 cubes of edge-length 1/6, for a total 
volume of 2/27, which, as expected, is 4/9 of vo^S'i) = 1/6. 

Discretizing this new shape with scale factors .si < S2 < ■§3 gives a subset of of 
size I6S1S2S3 which gives a lattice tiling of Z"^; this subset is included in 5"^, where k = 
si + 2s2 + 3s3 — 3. Another simple optimization shows that, for any given /c > 3, we get a 
subset of S'f^ which tiles and has size 16030406, where Oj = [{k + i)/6\. For large k (in fact, 
for all k > 30) , this new lattice gives a better covering of Z^ by S'j^ than the preceding one 
did, but for smaller k the preceding one sometimes does better. 

Aguilo, Fiol, and Garcia [1] also work with this shape, but discretize it in a rotationally 
symmetric way rather than in each dimension separately; the Cayley graphs they obtain 
are slightly larger than the graphs of size 16030405 given above, but still of the form 
(2/27)A;3 + 0(A;2). 

In order to see whether these lattice coverings give close-to-optimal Cayley graphs, the 
authors performed a computer search for the best (smallest-diameter) directed Abelian 
Cayley graphs on three generators. This extends similar computations performed by 
Aguilo, Fiol, and Garcia [1] and by Fiduccia, Forcade, and Zito [14]. The latter paper 
also contains a useful upper bound: an Abelian Cayley digraph on three generators with 
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diameter k has size at most 3{k + 3)^/25. This improves the obvious upper bound l^fel 
when k > 7. 

Comparing the above figures with the output from the authors' computations gives a 
shght surprise: the best cychc groups do substantially better than the groups from the 
above coverings. The data are shown in Table 2; here 'FFZ' is the Fiduccia-Forcade-Zito 
upper bound, 'Impr.' refers to the larger of the sizes obtained from the two improved 
constructions above, 'AFG' is the size attained by the Aguilo-Fiol-Garcia construction, 
and the remaining columns are analogous to those of Table 1. The computations were run 
on Abelian groups of sizes up to and including 4871; this means that the entries marked 
with an asterisk in the column (for which the FFZ bound is greater than 4871) have 
not been completely proven optimal, but it is extremely likely that they are. 

Note that in three cases, k = 7,31,33, the best cyclic Cayley graph was not achieved 
using 1 as one of the generators. If one is required to use 1 as a generator (which may 
be useful when actually building the corresponding loop network), then the best one can 
do is size 78 for A; = 7 (with generators 1,6,49), size 3178 for A; = 31 (with generators 
1, 386, 1295), and size 3794 for A; = 33 (with generators 1,469, 2094). 

There is one other difference between this case and the undirected case: here there 
are values of k for which one can do better with general Abelian groups than with cyclic 
groups. The improved values obtained from non-cyclic groups are shown in Table 3. A 
number of these optimal graphs are actually obtained by applying Proposition 5(b) to 
smaller Cayley graphs; for instance, the Abelian graph for A; = 17 is obtained this way 
from the cyclic graph for k = 7, which is the reason that these two graphs give exactly the 
same real-covering efficiency (.504). 

The values in the n'^ column of Table 2 are so much larger than those in the preceding 
two columns that it is clear that the real lattices used for the preceding columns were not 
optimal. This is made explicit in the last column of the table, which gives the efficiency 
of the real lattice covering obtained from the computed integer lattice covering via Propo- 
sition 4(b). For A; = 1 and A; = 3 these coverings are just (scaled versions of) the two 
coverings we explicitly constructed above; but later coverings obviously do substantially 
better. 

The best real covering obtained from these computations is that for k = 7, with effi- 
ciency .504. As in the undirected case, we can reconstruct generators for the lattice from 
the given generating set 2,9,35 for Zg^; after simplification, the resulting generating vec- 
tors are (—2, 2, 2), (3, —3, 3), and (4, 3, —1). We now have a computer-assisted proof that 
the lattice generated by these vectors gives a lattice covering of R"^ by 5*^0! but one can 
obtain useful extra information (as well as, perhaps, more satisfaction) by proving this 
directly. 

Proposition 17. Let Ly be the lattice in generated by the vectors (—2, 2, 2), (3, —3, 3), 
and (4, 3, -1); then S[o + L'^ = R^. 
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k 


\S'k\ 


FFZ 


Toroidal 


Impr. 


AFG 


< 


Generators 


/ / 1 /"If 1 

n'c/lS'kl 


</vol(6'l,+3) 





1 




1 




1 


1 




1 


.222222 


1 


4 




2 


4 


4 


4 


1,2,3 


1 


.375000 


2 


10 




4 


4 


7 


9 


1,3,4 


.900000 


.432000 


3 


20 




8 


16 


16 


16 


1,4,5 


.800000 


.444444 


4 


35 




12 


16 


19 


27 


1,4, 17 


.771428 


.472303 


5 


56 




18 


32 


31 


40 


1,6,15 


.714286 


.468750 


6 


84 




27 


32 


50 


57 


1,13,33 


.678571 


.469136 


7 


120 


120 


36 


48 


56 


84 


2,9,35 


.700000 


.504000 


8 


165 


159 


48 


72 


86 


111 


1,31,69 


.672727 


.500376 


9 


220 


207 


64 


128 


128 


138 


1,11,78 


.627273 


.479167 


10 


286 


263 


80 


128 


134 


176 


1,17, 56 


.615385 


.480655 


11 


364 


329 


100 


144 


182 


217 


1,13,119 


.596154 


.474490 


12 


455 


405 


125 


192 


243 


273 


1,14,153 


.600000 


.485333 


13 


560 


491 


150 


256 


252 


340 


1,90,191 


.607143 


.498047 


14 


680 


589 


180 


288 


333 


395 


1,35,271 


.580882 


.482394 


15 


816 


699 


216 


432 


432 


462 


1,29,97 


.566176 


.475309 


16 


969 


823 


252 


432 


441 


560 


1,215,326 


.577915 


.489867 


17 


1140 


960 


294 


500 


549 


648 


1,76,237 


.568421 


.486000 


18 


1330 


1111 


343 


576 


676 


748 


1,41,147 


.562406 


.484613 


19 


1540 


1277 


392 


600 


688 


861 


1,27,463 


.559091 


.485162 


20 


1771 


1460 


448 


768 


844 


979 


1,22,351 


.552795 


.482781 


21 


2024 


1658 


512 


1024 


1024 


1140 


1,45,196 


.563241 


.494792 


22 


2300 


1875 


576 


1024 


1036 


1305 


1,246,1030 


.567391 


.501120 


23 


2600 


2109 


648 


1024 


1228 


1440 


1, 126,415 


.553846 


.491579 


24 


2925 


2361 


729 


1280 


1445 


1616 


1,56,257 


.552479 


.492608 


25 


3276 


2634 


810 


1372 


1460 


1788 


1,154,1452 


.545788 


.488703 


26 


3654 


2926 


900 


1600 


1715 


1963 


1,90,780 


.537219 


.482923 


27 


4060 


3240 


1000 


2000 


2000 


2224 


1,425,704 


.547783 


.494222 


28 


4495 


3574 


1100 


2000 


2015 


2442 


1,964,1372 


.543270 


.491826 


29 


4960 


3932 


1210 


2048 


2315 


2693 


1,39,942 


.542944 


.493103 


30 


5456 


4312 


1331 


2400 


2646 


2920 


1,540,831 


.535191 


.487520 


31 


5984 


4716 


1452 


2400 


2664 


3220 


7, 30, 2277 


.538102 


.491553 


32 


6545 


5145 


1584 


2880 


3042 


3591* 


1,1519,2031 


.548663 


.502531 


33 


7140 


5598 


1728 


3456 


3456 


3850* 


2,475,1177 


.539216 


.495113 


34 


7770 


6078 


1872 


3456 


3474 


4191* 


1,748,2652 


.539382 


.496437 


35 


8436 


6584 


2028 


3456 


3906 


4468* 


1,353,2789 


.529635 


.488555 



Table 2. Best directed Cayley graphs of cyclic groups, three generators. 
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7 

k 




Group 


Generators 




</vol(5^+3) 


12 


279 


Z93 X Z3 


(1,0), (9,1), (10, 2) 


.613187 


.496000 


17 


672 


Zigg X Z2 X Z2 


(2, 1,0), (9, 0,0), (35, 0,1) 


.589474 


.504000 


18 


752 


Z188 X Z4 


(1,0), (13, 2), (14,1) 


.565414 


.487204 


19 


888 


Z222 X Z2 X Z2 


(1,0,0), (142, 1,0), (180, 0,1) 


.576623 


.500376 


26 


1980 


Z330 X Zg 


(1,0), (123, 2), (234, 3) 


.541872 


.487105 


27 


2268 


Z252 X Z3 X Z3 


(2, 0,0), (9, 1,0), (35, 0,1) 


.558621 


.504000 


28 


2448 


Zsie X Z3 


(1,0), (427, 0), (564,1) 


.544605 


.493035 


29 


2720 


Zeso X Z2 X Z2 


(1,0,0), (191, 1,0), (90, 0,1) 


.548387 


.498047 


30 


2997 


Z333 X Z3 X Z3 


(1,0,0), (31, 1,0), (180, 0,1) 


.549304 


.500376 


35 


1500* 


Z:5()() X Zi,-j 


(1.0). (3.1). (21 1.7) 


.533128 


.192051 



Table 3. Best directed Cayley graphs of Abelian groups, three generators. 
Proof. First, note that the following vectors are in L^: 



Vi = 


(-2,2,2) 




Vs = 


(1,-1,5) 


= Vi + V2 


V2 = 


(3,-3,3) 




V9 = 


(-1,1,7) 


= 2vi + V2 


V3 


(4,3,-1) 




VlO 


(3,4,-6) 


= V3 - Vi - V2 


V I = 


(6,1,-3) 


= V3 - Vi 


V| 1 = 


(-7,7,1) 


= 2vi - V2 


V5 = 


(5,-5,1) 


= V2 - Vi 


V12 = 


(-5,-2,8) 


= 2vi + V2 - V3 


V6 = 


(1,6,-4) 


= V3 - V2 


Vi3 = 


(-1,8,-2) 


= Vi - V2 + V3 


V7 = 


(2,5,1) 


= Vi + V3 


Vi4 = 


(8,-1,-5) 


= V3 - 2vi 



For each vector = (r, s, t), we have l<r + s + t<8; hence, the translated tetrahedron 
S'lQ + Vj intersects the plane x + y + z = 10. In fact, the intersection is a triangle whose 
vertices have coordinates (10— s— s, t), (r, 10— r— t,t), and (r, s, 10— r— s). 

Figure 4 shows the upper face of S[q. For each i < 14, it indicates which part of this 
face is covered by the translate S[q + v^. (The face is divided up into unit triangles, each 
of which is labeled by the value(s) of i for which S[q + covers that triangle.) Clearly 
each unit triangle is labeled, so the translates S[q + Vj cover the entire upper face of S[q. 

In fact, since each Vj has coordinates summing to at least 1, the smaller translates 
Sq + i < 14, cover the upper face of 5'(q. For any x in this upper face, there is an i 
such that X e S'g + Vi, so S'i + X C 5( + + Vj = S[q + Vj. Since S'n is the union of 
S[q and the sets + x for x in the upper face of S'-^q, we have S[-^ C S[q + L^, and hence 

We now prove by induction that, for all integers A; > 10, 5^ + L'^ C S[q + L7. The case 
A; = 10 is trivial. If it is true for /s, then 

"^fe+i + -^7 — "^1 + "^fe + -^7 ^'5^1 + S'lQ + = S'li + C S[q + L7, 
so it is true for A; + 1. 
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(0,0,10) 




Figure 4. Coverage of one face of S[q under L^. 

Finally, for any y G R^, there is a member w of L'^ such that the coordinates of y — w 
are all positive (e.g., let w be a large multiple of — V7). Then y — w e S'^ for some k, so 
y e S'f^ + L'^ and hence y G + ^7- Therefore, S[q + L'^ = 11^. ■ 

This covering has efficiency det(vi, V2, V3)/ vol(S'(o) = -504. Hence, Corollary 10(b) 
gives: 

Corollary 18. For all k, there is a directed Cayley graph of an Abelian group on three 
generators which has diameter k and size at least O.Q84k^ + 0{k^). ■ 

We can now use the method of (the real version of) Lemma 3 to get a fundamen- 
tal region T7 C S[q for the lattice Ly. (Recall that any such region must have volume 
det(vi, V2, V3) = 84.) To do this, just start with S'iq, look at each of the vectors Vj 
{i < 14) defined above, and delete those points of S'iq which lie in S[q + v^. (We have 
Vi y for all i if ~< orders vectors primarily by the sum of coefficients.) What is left is 
the set shown in Figure 5; since this is the union of 84 unit cubes, we know that there is 
no need to subtract further translates S[q + w. This set was obtained independently by 
Fiduccia, Forcade, and Zito [14]. 
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Figure 5. A subset of the tetrahedron S'iq which tiles space. 

So this set Ty gives a lattice tiling of using Ly. This tiling is quite unusual; the 
translates of fit together in a peculiar way, seeming to wind around each other. One 
interesting fact is that each translate of T7 is adjacent to (i.e., shares a boundary segment 
of positive area with) 28 other translates, a surprisingly high number. (Ty itself is adjacent 
to + Vj and - Vj for i < 14.) 

In many of the tilings we constructed explicitly, there was a polycube fundamental region 
like T7, but there was also an alternate fundamental region which was convex; for instance, 
for the optimal covering of by right triangles, one could use either an L-tromino or a 
hexagon as the fundamental region. Clearly Ly has convex fundamental regions (e.g., its 
Voronoi regions), but it turns out that they are unsuitable for the current problem: 

Proposition 19. There is no convex fundamental region for the lattice L'j included within 
the tetrahedron S[q. 

Proof. Since Ty gives a lattice tiling of R^ by Ly, every point of R"^, except for those lying 
on boundaries of the tiling, can be translated by a vector in Ly to a unique point of Ty. 
In particular, if we look at the part of S[q lying outside Ty, then we can break it up into 
finitely many parts (in fact, just cut it along the integer translates of the three coordinate 
planes) which can be translated in a unique way by members of Ly so as to lie within Ty. 

If one does this, one finds that there are parts of Ty which do not get covered by 
translates of parts of S[q \ Ty. Most of these uncovered parts look like inverted copies 
of S'l (i.e., translates of — /§(), although there are some larger ones. In particular, the sets 
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(1, 1, 1)—S'i, (8, 1, 1)—S'i, (1, 8, 1) — S[, and (1, 1, S) — S[ are not covered by such translates. 
This imphes that each of those four sets is disjoint (except for boundaries) from all of the 
translates 5^Jo + w for w G Ly\{0}. It follows that any fundamental region for included 
within S'lQ must include all four of these sets. 

If the fundamental region is also convex, then it must contain any convex combinations of 
points in those four sets; in particular, it must include the sets (3, 1, 1)— and (1, 3, 3)—S'i. 
But (1, 3, 3) = (3, 1, 1) + Vi, so the vi-translate of the region overlaps the region itself in a 
set of positive volume, which is impossible for a fundamental region of Ly. Therefore, no 
fundamental region of within S[q can be convex. ■ 

There is no obvious reason why the lattice Lj should be exactly optimal for a lattice 
covering of by the tetrahedron S[q. In the case of undirected graphs on three generators, 
the lattices obtained for each k were not optimal, but were closer and closer approximations 
to the lattice Lbcc, which does appear to be optimal. One would expect something similar 
to occur in the directed case, but it does not; the real lattice efficiencies in the last columns 
of Tables 2 and 3 go up and down irregularly and do not (so far) exceed the value .504 
attained by L'-^. 

Given this, it seems reasonable to examine and try to adjust it slightly in order 
to improve its efficiency; there should be some locally optimal lattice which Ly is an 
approximation to, and we would like to find it. Quite surprisingly, it turns out that no 
adjustment is necessary. Just before submitting the present paper, the authors found a 
recent paper of Forcade and Lamoreaux [15] proving this same result by a method slightly 
different from that presented here. 

Theorem 20 (Forcade and Lamoreaux [15]). Among those lattices L for which S[q + L = 
R^, the lattice Lij is locally optimal. 

Proof. We use the same methods as for Theorem 16. Recall the vectors vi, . . . , V14 from 

Proposition 18. The vectors vi, V2, V3 generate Ly; a nearby lattice L will be generated by 
nearby vectors = (ai,6i,ci), V2 = («27^25C2), and Vg = (03,63,03). Again concatenate 
v[, V2, V3 and vi, V2, V3 to get v' and v in R^. Let i^(v') be the determinant of the matrix 
with rows V2, V3; then we have F{v) — 84, and we want to see that -F(v') < 84 for any 
other v' near v for which the corresponding lattice L satisfies S[q + L = R^. We compute 
that the gradient of F at the point v is g = (-6, 15, 21, 8, -6, 14, 12, 12, 0). 

Referring back to Figure 4, we see that the point (1,1,8) is on the boundary of S[q + Vj 
for i = 8,9, 12, as well as on the boundary of .S^o itself; one can check that no other 
Ly-translate of S[q is near this point. The nearby lattice L contains points 0, Vg = 
v[ + V2, Vg = 2v[ + V2, and v'12 = 2vi + V2 — V3. For any small positive s, the point 
(oi + a2 — 2bi + 62 — £, 2ci + C2 — cs — e), which is near (1, 1, 8), will not be in Vg + S[q 
because its first coordinate is smaller that that of Vg. By looking at second and third 
coordinates respectively, we see that this point is not in Vg + S'iq or v'12 + ^'10 either. 
Hence, in order to have S[q + L = R^, this point must be in S[q itself, so we must have 



ai + aa + 2bi + 62 + 2ci + C2 - C3 - 3£ < 10. 
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Since e can be arbitrarily small, we need 

ai + a2 + 26i + 62 + 2ci + C2 - C3 < 10 
in order to have S'iq + L = R^. So we have the constraint 

ui • v' < 10, where Ui = (1, 2, 2, 1, 1, 1, 0, 0, -1). 

The same reasoning applied at the points (1,2,7), (3,2,5), (5,2,3), (5,3,2), (4,4,2), 
(3, 5, 2), (2, 6, 2), (1, 7, 2), (8, 1, 1), (6, 3, 1), (3, 6, 1), and (1, 8, 1) gives the constraints 



U2 


v' 


< 10, 


where 


U2 = 


(1,1,2,1,0,1,0,0,0), 


U3 


v' 


< 10, 


where 


Us = 


(0,1,1,1,0,1,0,0,0), 


U4 


v' 


< 10, 


where 


U4 = 


(-1,1,0,1,0,1,0,0,0), 


U5 


v' 


< 10, 


where 


U5 = 


(-1,0,1,1,0,0,0,1,0), 


U6 


v' 


< 10, 


where 


U6 = 


(0,-1,1,0,-1,0,1,1,0), 


U7 


v' 


< 10, 


where 


U7 = 


(-1,1,1,-1,0,0,1,1,0), 


Us 


/ 

V 


< 10, 


where 


Us = 


(1,0,1,0,-1,0,1,1,0), 


ug 


/ 

V 


< 10, 


where 


Ug = 


(0,2,1,-1,-1,0,1,0,0), 


UlO 


v' 


< 10, 


where 


UlO = 


(-2,-1,-1,0,0,1,1,1,0), 


Ull 


v' 


< 10, 


where 


Ull = 


(-1,0,-1,0,0,1,1,1,0), 


U12 


v' 


< 10, 


where 


U12 = 


(-1,0,1,-1,-1,0,1,1,1), 


Ul3 


v' 


< 10, 


where 


Ul3 = 


(0,1,2,-1,-1,-1,1,1,0). 



Again note that all thirteen of these inequalities are satisfied with equality when v' = v. 
Hence, they can be rewritten as Uj • (v' — v) < for z = 1, 2, . . . , 13. 

One can easily check that the vectors ui, U2, U4, U5, ug, uy, Us, Uio are linearly indepen- 
dent; their common null space is generated by the vector w = (0, 0, 0, 1, 1, —1, 2, —1, 1). 
(The other five vectors are also orthogonal to w, so they are linear combinations of the 
eight listed above.) Also, we have 

g = Ui + 4.8U2 + 6.4U4 + U5 + 1.6U6 + 3.2U7 + 3.4U8 + Ug + l.Suio + U12. 

Let C be the closed cone consisting of all vectors t in the subspace spanned by Ui, . . . , U13 
such that Uj ■ t < for all i < 13. Then the above equations imply that g ■ t < for all 
t in C, and equality can hold only when t = 0. Hence, as in Theorem 16, there is a 
neighborhood t/ of g such that, for any g' in U and any nonzero t in C, g' ■ t < 0. 

We can compute that, for any real number r, the determinant for the lattice given by 
V + rw is F(v + rw) = 84 — 12r^. Clearly this is at most 84, with equality only when 
r = 0. 
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Now, any vector v' close to v can be expressed as v + ti +t2 where ti is a small multiple 
of w and t2 is a small linear combination of the vectors Ui, . . . , U13. The reasoning from 
Theorem 16 shows that t2 must be in C if S'^q + L = R^. Also as in that Theorem, we find 
that F(v + ti) < F(v) with equality only when ti = 0, and F(v + ti + < F{v + ti) 
with equality only when t2 = 0. Therefore, i^(v') < i^(v), with equality holding only 
when v' = v. So v gives a local maximum of F, as desired. ■ 

This and the computational evidence make it plausible that Lj actually gives an optimal 
lattice covering of R"^ by ^^q, and hence that the asymptotic formula in Corollary 18 is 
optimal. 

More than three generators. In higher dimensions, analogues of many of the preceding 
constructions exist, but they do not produce lattice coverings as efficient as one would hope 
for. 

For lattice coverings with the d-dimensional dual cube, one can use the d-dimcnsional 
body-centered cubic lattice (the set of vectors in whose coordinates are all odd or all 
even). By the same argument as for the three-dimensional case, this lattice gives a lattice 
covering of R'' by 8^/2- The efficiency of this covering is 2^~^ / vol{Sd/2) ~ 2'^~^d\/cfi, 
which is 2*^"^ times the efficiency of the covering using the ordinary cubic lattice Z^. 

As usual, the Cayley graph corresponding to this lattice is a twisted toroidal mesh. For 
a given number m, one can connect the elements of 'i'^rn ^ '^rn as in an ordinary toroidal 
mesh, except that the wraparound connections for the last coordinate are twisted along 
all of the other coordinates: (xi, . . . , x^-i, to— 1) is connected to (xiiTO, . . . , x^-iiTO, 0). 
This gives a graph of diameter [dm/2\ and size 2^~^m^, which is about 2^~^ times as 
large as the best ordinary toroidal mesh of this diameter. 

One can optimize this slightly. Given the dimension d and the desired diameter k, let 
q and r be the quotient and remainder when 2A; + 1 is divided by d; wc assume k is large 
enough that q > 0. Then a good lattice L to use is the body-centered cubic lattice above 
scaled up by a factor q + 1 in each of the first r coordinates and a factor q in the remaining 
d—r coordinates. The resulting Z'^/L is isomorphic to (Z2q+2 ^ '^2q^)/H with the canonical 
generators, where H is the two-element subgroup {0, (g+1, . . . , g+1, g, . . . , q)} (there are 
r g+l's); it can be laid out as a twisted toroidal mesh on Z2q_|_2 x ^2"^"^ x Zg or on 
^2g+2 ^ '^2q^ ^ Zg+i- If 9 is even and r > 0, this Cayley graph is isomorphic to that of 
Zg_|_i X Z2~^2 ^ '^tq'^ with the generators e2, . . . , e^; and (1, . . . , 1, g, . . . , g) with r I's; if 
g is odd or r = 0, then it is isomorphic to the Cayley graph of Z2g_j.2 ^ '^2q^~^ ^ ^9 with 
generators ei, . . . , e^-i and (g+1, . . . , g+1, 1, . . . , 1) with r g+l's. The size of this graph 
is slightly larger than the size of the cyclic Cayley graph constructed by Chen and Jia [5] , 
but the ratio of the two sizes tends to 1 for large k. 

For the directed case, we must consider lattice coverings by d-simplices; as usual, by 
affine invariance, it doesn't matter which simplex is used. One can show that a lattice for 
covering with a given d-simplex is given by the following generating vectors: for each face 
of the simplex, take a vector which is twice the vector from the centroid of the simplex to 
the centroid of that face. (This gives d+1 vectors, but they sum to 0, so just take d of 
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them.) The efficiency of this covering works out to be d\2'^ / {d'^{d + 1)). 

Unfortunately, in both cases, the efficiency decreases exponentially with d: by Stirling's 
formula. 



This seems to be the case for all known explicitly constructed lattice coverings by these 
shapes (and by spheres). 

On the other hand, in 1959 Rogers [18] gave a nonconstructive proof that there exist 
much more efficient lattice coverings by these shapes (or by any convex body) in high 
dimensions; and he gave an even better result for the case of spheres. More recently 
Gritzmann [16] extended the latter result to apply to any convex body with a sufficient 
number of mutually orthogonal hyperplanes of symmetry. (The number required is quite 
small: only |_log2ln(iJ +5.) Gritzmann's result states that there is a constant c (not 
depending on d or on the convex body) such that, for any convex body K in R*^ with the 
above number of mutually orthogonal planes of symmetry, there is a lattice covering of R'' 
by K with density at most cd{lii d)^'^^°^^ ^ . 

The regular dual d-cube and the regular d-simplex do have the required symmetry 
planes for large enough d. This is clear for the dual d-cube; it has the same d orthogonal 
planes of symmetry as the d-cube it is dual to. For the regular cZ-simplex, note that the 
perpendicular bisector of an edge is a hyperplane of symmetry, and that edges which do 
not share a vertex have orthogonal directions (the easiest way to see this is to look at the 
regular d-simplex in R*^"*"^ whose vertices are ei, . . . , ed+i, and take dot products), so one 
can find \d/2~\ mutually orthogonal hyperplanes of symmetry. Therefore, we get lattice 
coverings of the specified density for large d, and by adjusting the constant c we can make 
the bound apply for all d (for these two particular shapes). Therefore, letting c = c~^, we 
can use Corollary 10 to get: 

Theorem 21. There is a constant c > (not depending on d or k) such that, for any fixed 
d> 1 and for all k, there exist undirected Cayley graphs of Abelian groups on d generators 
having diameter < k and size at least 




and 




-k' + Oik"-'), 



d!d(lnd)i+iog: 



and there exist directed Cayley graphs of Abelian groups on d generators having diameter 
< k and size at least 
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The coverings produced by this method are probably fairly strange. We seem to have 
run into this already in three dimensions, for the directed case; for the undirected case it 
apparently happens later. 

Layouts with short wires. The obvious way to lay out a toroidal mesh is as a rectangular 

array with mesh connections between adjacent nodes in the array and with long wires 
connecting opposite ends of the array; these long wires may cause communications delays. 
However, there is a standard trick for rearranging the layout so as to remove the need 
for long wires. In the one-dimensional case, instead of placing the nodes in the order 
1,2,3, ... ,n (where i is connected to z + 1 for i = 1,2, . . . , n—1 and n is connected to 1), 
one can place them in the order 1, n, 2, n—1, 3, n—2, . . . ; then the maximum required wire 
length is only twice the mesh spacing. In higher dimensions, one can apply the same trick 
to each dimension separately, and again the required wire length is twice the mesh spacing. 

It is not immediately obvious that this interleaving trick can be applied to twisted 
toroidal meshes; a simple interleaving in each dimension would not make the twisted 
cross-connections short. But it is possible to get short-wire layouts for the twisted meshes 
in a similar way. One approach is to perform the interleaving twice on the long dimensions 
of the mesh; for instance, if the mesh has length 16 in one of the long dimensions, then 
the nodes would be arranged in the order 

1, 9, 16, 8, 2, 10, 15, 7, 3, 11, 14, 6, 4, 12, 13, 5. 

Then wires in this dimension would have length at most 4 times the mesh spacing. Now, 
when one does a single interleaving on the short dimension, the twisted cross-connections 
become short as well. 
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Figure 6. Short-wire layout for a twisted toroidal mesh. 



Another method is shown in Figure 6. Here the idea is to modify the original arrange- 
ment (a) by shearing the mesh (rotating the z'th level in the short dimension by z — 1 units 
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in each of the long dimensions), as shown in (b), so that the twisted cross-connections 
become (almost) straight. Then one can do an ordinary interleaving in each dimension to 
get the result shown in (c). This gives a layout in which the maximum wire length is 2\fd 
times the mesh spacing. 

Some of the other Abelian Cayley graphs we have considered can be treated similarly, 
especially the ones which differ only slightly from twisted toroidal meshes. For the optimal 
two-generator undirected Abelian Cayley graph, if one starts with the almost-rectangular 
layout shown in Figure 1 (a (fc-l-l) x (/c-|-l) square next to a, k x k square), and performs 
a shear as in Figure 6, then the result is a 2A; x (k+l) rectangle with one node left over; 
this can then be interleaved to get a short- wire layout. A more difficult case is the two- 
generator directed graph from Theorem 13 and Corollary 14; here one can start with the 
natural L-shaped layout and perform shears on separate parts to obtain a rectangle with 
dimensions (a + 2b) x a (made up of three subrectangles with different shear patterns) 
where the necessary cross-connections are almost straight across, and hence interleaving 
will give a good layout. 

Generators of order 2. The undirected Cayley graphs produced so far all have even 
degree (twice the number of generators). If one is interested in undirected Cayley graphs 
of odd degree, one will have to use a generator of order 2. 

Using d unrestricted generators plus one order-2 generator, one can get an undirected 
Abelian Cayley graph of a given diameter which is about twice as large as one can get using 
d unrestricted generators alone. More precisely, if na{d, k) is the size of the largest Abelian 
Cayley graph of diameter k using d generators and n+((i, k) is the size of the largest such 
graph using d generators plus one order-2 generator, then 

2n„(d, k-1) < n+{d, k) < 2n„(d, k). 

To see this, first let G be generated by gi, . . . ,g(i and p where p has order 2. If G has 
diameter at most k using these generators, and H is the subgroup of size 2 generated by p, 
then G/H is generated by the images Qi + H for 1 < i < d with diameter at most k, so 
\G/H\ < na{d,k), so \G\ < 2na{d,k); hence, n'^{d,k) < 2na{d,k). On the other hand, if 
G is generated by (71, . . . , (/^ with diameter at most k — 1, then G x Z2 is generated by 
{gi, 0) for i — 1, . . . ,d and (0, 1), and has diameter at most k using these generators; this 
shows that 2na{d, k — 1) <n'^{d,k). 

We can also study ?i+(d, k) using the same methods that were used for na{d,k). The 
appropriate free (universal) group to use here is the group Z*^ x Z2, with the canonical 
generators (e^, 0) for i = 1, . . . ,d and (0, 1). The set of elements of this group which can be 
written as a word of length at most k in the generators is precisely Wk — {SkXO)U{Sk-i x 1). 
For any Abelian group G with generators gi, . . . ,gci and p (p of order 2), there is a unique 
homomorphism from x Z2 to G taking (cj, 0) to gi and (0, 1) to p; the Cayley graph of G 
using these generators has diameter at most k if and only if the homomorphism maps 
onto G. So the obvious upper limit for the size of G is \Wk\ = \Sk \ + |5'jt-i|. 

We are now led to study quotient groups (Z*^ x Z2)/N where is a (normal) subgroup 
of X Z2 of finite index; we want such an N of largest possible index such that Wk + N = 
Z'^ X Z2. 
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One simple possibility is that C x {0}; in this case the resulting group is just 
(Z'^/Nq) X Z2 where N = Nq x {0}. It is easy to see that the diameter of this group is 
precisely one more than the diameter of Z'^/Nq using the canonical d generators. 

Note that A'^o is a rf-dimensional lattice; let vi , . . . , be a list of generators for this 
lattice. Now let N' be the subgroup of Z'^ x Z2 generated by (v^, 1) for i = 1, . . . , d. Then 
we have 

IZ"^ X Z2 : A^'l = 2 \Z'^ : A^ol = IZ'^ x Z2 : iV| . 

Furthermore, the diameter of (Z'^ x Z2)/A^' is at most one more than the diameter of Z'^/Nq, 
which means that it is no larger than the diameter of (Z*^ x Z2)/N. 

This shows that, when trying to determine k), we may restrict ourselves to study- 

ing subgroups AT of Z** x Z2 of finite index which are not included in Z'^ x {0}. 

So choose g e Z*^ such that (g, 1) e N. The subgroup ATn (Z*^ x {0}) is of index 2 
in N and hence of finite index in Z'^ X Zo. So we have iV n (Z<^ X {0}) = L x {0} for some 
(i-dimensional lattice L. Note that (2g, 0) = 2(g, 1) G N, so 2g G L. (Normally g will not 
be in L; if g G L, then (g, 0) G A^, so (0, 1) G A^, so the order-2 generator collapses to the 
identity in the quotient group.) Also, we have [Z*^ x Z2 : A?"! = \Z'^ : L\. 

It is now easy to see that 

{Wk + N)n {Z^ X {0}) = {{Sk + L)U {Sk-1 + g + L)) X {0}. 

Hence, in order to have Wk + A'" = Z'^ x Z2 , it is necessary to have 

(*) {Sk + L)U{Sk-i+g + L)^Z'^. 

This necessary condition is also sufficient, because 

(Wk + N)n (Z^ X {!}) = {{Sk-1 +L)u{Sk + S + L)) x {1} 

= {{{Sk + L) U {Sk-1 + g + L)) + g) X {!}. 

So our goal is to find such a lattice L and extra generator g (with 2g G L) so that (*) is 
satisfied and jZ*^ : L| is as large as possible. 

We are now ready to consider specific values of d. As usual, the case d — 1 is easy. The 
maximal possible value of \Z : L\ is \Wk\ — 4/c, and this value is attained by letting L be 
generated by the element 4/c, with g = 2k. This leads to the cyclic Cayley graph on the 
group Z^k with unrestricted generator 1 and order-2 generator 2k. 

For d = 2 we have a situation very similar to that in Figure 1 (lattice coverings with 
Aztec diamonds), but not identical because we must use two different shapes. The upper 
bound on |Z^ : L\ is \Wk\ — ^k"^ + 2. For k — 1 this bound is actually attainable; it leads 
to the Cayley graph from Zq with unrestricted generators 1 and 2 and order-2 generator 3. 
But for A; > 1 the pieces Sk and Sk-i do not fit together well enough to give a perfect 
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Figure 7. Lattice covering of T? using S}. and S\~-\ (shown for A; = 3). 

tiling. The best one can do is the lattice L generated by (2A; + 1, 1) and (—1, 2A; — 1), with 
the extra generator g = (A;, A;), as shown in Figure 7. 

The graph of diameter k resulting from this covering is the Cayley graph of the cyclic 
group Z4fc2 with unrestricted generators 1 and 2 A; — 1 and order-2 generator 2A;^. One can 
get another Cayley graph of this size by using the lattice generated by (2A;, 0) and (0, 2A;), 
but this graph will not be cyclic; it comes from the group T^ih x '^ik- 

The outlined shape in Figure 7 (a (2A; + 1) x (2A; — 1) rectangle with one extra point) 
is a fundamental region which is convenient for an actual layout of nodes in a network. 
Note that a 2/c x 2/c rectangle (or for that matter, any rectangle with both sides greater 
than 1) will not work as a fundamental region for this lattice. The alternative lattice in 
the previous paragraph does allow one to use a layout which is a 2A; x 2A; rectangle; in 
fact, this is just a toroidal mesh. However, in either case one will have to make the extra 
connections specified by the order-2 generator. 

When one moves to (i = 3, it becomes harder to get optimal results, so again the 
authors resorted to a computational search. For k = 1 the best graph is the Cayley graph 
of Zg with unrestricted generators 1,2,3 and order-2 generator 4; for A; = 2 the best is 
Z26 with unrestricted generators 1,2,8 and order-2 generator 13. For 3 < A; < 10 the 
optimal results, like those for three generators alone, form a pattern of period 3, as shown 



THE DEGREE-DIAMETER PROBLEM FOR ABELIAN CAYLEY GRAPHS 



41 



in Table 4. (Again the best graphs are aU cychc. This time the parameter a is defined to 
be the integer nearest 2k/3.) 



k mod 3 





1 


2 


a 


2A;/3 


(2/e+l)/3 


(2A;- l)/3 


Lattice 
f?pn era.tors 


(2a, 1,-1) 
(-l,2a,-l) 
(1,1,2a) 


(2a -1,-1,0) 

(l,2a,-l) 
(0,1,2a- 1) 


(2a + 1,-1,0) 
(l,2a,-l) 
(0,l,2a+l) 
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I /~i /~i /~i 1 1 

l^ci, a, a — i-j 


^^a ± , a, aj 


Cychc 
group size 


Uk^ + 108A; 
27 


64A;3 + 60A; - 16 
27 


64A;3 + eO/c + 16 
27 


Unrestricted 
generators 


1 

Aa^ -2a^ + 2a-l 
4a^ - 2a^ + 4a - 1 


1 

2a- 1 
4a^ - 2a + 1 


1 

2a + 1 
4a^ + 2a + 1 


Order-2 
generator 


4a3 + 3a 


4a3 - 4a2 + 3a - 1 


Aa^ + 4a^ + 3a + 1 



Table 4. Best undirected Cayley graphs of cyclic groups of diameter k using 
three generators plus one order-2 generator (3 < /c < 10). 

We can now apply the methods in the proof of Theorem 15 to show: 

Theorem 22. For each k > ?>, the cyclic undirected Cayley graph using the group and 
generators specified in Table 4 has diameter k. ■ 

The authors again conjecture that these are actually the optimal such Abelian Cayley 
graphs for all A; > 3, not just for 3 < A; < 10. 

For d > 3 one can get reasonably good results by letting L be approximately a cubic 
lattice, with g near the center of one of the cubes; this makes L U (g + L) approximately 
a body-centered cubic lattice. Again, though, the efficiency of this covering decreases 
exponentially with d; one can do much better using the results of Gritzmann [16]. 

One can also consider the possibility of using more than one generator of order 2. For 
instance, one could look at Cayley graphs of degree 2d + 2 obtained by using d unrestricted 
generators and two generators of order 2. 

However, this is not going to he helpful if one wants to construct large undirected 
Cayley graphs of a given degree and diameter, at least if the diameter is substantially 
larger than the degree. For instance, suppose that the degree is fixed as 2d + 2. If one 
uses d unrestricted generators and two order-2 generators, then the size of the resulting 
undirected Abelian Cayley graph of diameter k is at most 4 times the number of points in 
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the (i-dimensional shape Sk- (More precisely, by the methods used above for one order-2 
generator, one gets a hmit of \Sk\ + 2\Sk-i\ + \Sk-2\-) This hmit is 0{k^), which is less 
than the 0{k'^~^^) one gets by using d + 1 unrestricted generators. The same argument 
shows that using more than two order-2 generators cannot be useful for large k; one gets 
larger graphs by replacing two order-2 generators with one unrestricted generator. 

If one is interested in the small-diameter case, though (especially when the diameter is 
less than or equal to the degree), then order-2 generators must be considered. The most 
extreme version of this would be to make all of the generators have order 2. This then 
becomes precisely the covering radius problem for binary linear codes; see the surveys by 
Cohen et al. [8,9] for more on this problem. The resulting graphs would be hyper cubes 
with additional diagonal connections to reduce the diameter. 

Conclusions. We have shown that one can construct Cayley graphs of Abelian groups 
which have substantially more vertices than traditional toroidal meshes, but retain certain 
desirable features. In particular, routing on the twisted toroidal meshes is easily described 
in almost the same manner as on toroidal meshes, and the twisted toroidal meshes host the 
discrete non-periodic orthogonal grids used in numerical calculations in exactly the same 
way that toroidal meshes do. In addition, we have shown how our d-dimensional meshes can 
be constructed with physical wire lengths that remain constant with increasing diameter 
(and increasing number of vertices) just as the corresponding toroidal meshes can. We 
have given results which are provably optimal in 2 dimensions and probably optimal in 3 
dimensions — the physically interesting cases. 

In the sequel to this paper, we will show that our methods can be extended to cover 
certain types of nilpotent groups. These groups yield graphs with cardinalities which still 
increase polynomially with diameter for a given degree, but with an exponent which is 
larger than in the Abelian case. One class of groups for which we obtain optimal results 
includes the groups discussed in Draper and Faber [12]. In particular, we show that the 
particular groups analyzed in that paper are not optimal for large diameters. 
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